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Abstract 

These lectures review some of the progress made in the quantitative understanding of 
B decays. The emphasis here is on applications of QCD using perturbative and non- 
perturbative techniques. In some cases, however, phenomenological models must at 
present be invoked to make meaningful comparison with data. The resulting picture 
is consistent with the standard model (SM) and this agreement is quantified in terms of 
the branching ratios, mixing probabilities, and lifetimes which measure the charge current 
and effective flavour changing neutral current transitions involving B hadrons. This, in 
turn, enables a determination of five of the nine elements of the quark mixing matrix. 
We discuss several proposals on improving the precision on the parameters of this matrix 
in forthcoming experiments. Issues intimately related to the quark mixing matrix such 
as the profile of the unitarity triangle and CP-violating asymmetries in B decays are 
discussed. In particular, we emphasize the role of rare B decays and B° - B° mixings in 
testing the SM quantitatively and in searching for physics beyond the SM. 
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1 Introduction 



The principal interest in the studies of B decays in the context of the standard model 
(SM) 0] lies in that they provide valuable information on the weak rotation matrix - the 
Cabibbo-Kobayashi-Maskawa (CKM) matrix M, [J. In fact, B decays determine five of 
the nine CKM matrix elements: V c b, V u b, Vtd, V ts , and Vtb- The dominant decays of 6-quark 
stem from the direct bcW~ coupling; then there are decay modes which stem from the 
CKM-suppressed buW~ coupling. These two classes represent the direct charged cur- 
rent (CC) transitions. The electromagnetic penguins and particle-antiparticle mixing(s), 
representing the so-called flavour changing neutral current processes (FCNC) which have 
been observed in B decays, are induced as higher order effects (loops) as the SM does 
not allow direct couplings of the form bsX or bdX, where X = 7, Z, H°, qq or a gluon. 
The effective induced couplings in the SM are governed by the GIM mechanism @] and 
are dominated by the intermediate (virtual) top quark contribution - the quark with the 
largest Yukawa coupling - through the transitions b — *> tW — > s and b — > tW — > d. Their 
quantitative measurements therefore provide information about the properties of the top 
quark, such as its mass and its weak mixing matrix elements V t d, V ts , and V t \,- 

To extract the CKM matrix elements from the hadronic transitions, one needs to im- 
plement the QCD perturbative corrections and calculate the hadronic decay form factors 
and decay functions for the inclusive and exclusive decays, respectively. A lot of work 
has gone into calculating the perturbative QCD corrections in B decays and this will be 
discussed in some detail here. The aspects having to do with non-perturbative physics 
are not yet under quantitative control, though important advances have been made and 
partial answers are available. In principle, non-perturbative aspects in B decays can all 
be calculated in the Lattice-QCD framework. In practice, the impact of this technique 
is limited due to the inadequacy of the present computer technology which restricts di- 
rect computation of the B decay properties involving the 6-quark with a typical mass 
of 0(5 GeV). However, useful information on some form factors and coupling constants 
has been obtained by simulation of the charmed hadron systems and extrapolating to the 
6-quark mass, often also using the constraints from the limiting behaviour of QCD in the 
rriQ — > 00 limit. There exist other non-perturbative theoretical tools such as the heavy 
quark effective theory HQET, the QCD sum rules, and the good old potential models, 
which have been put to good use in the quantitative analyses of experimental results in 
B decays. We shall review here some representative applications of each of these meth- 
ods. They, in particular the HQET techniques, have enabled us to determine the two 
mentioned matrix elements V c b and V u b- 

The CKM matrix elements Vu', i = d,s,b are, in principle, also measurable in the 
production and decays of the top quark ||. We note that first measurements of \V t b\ have 
been reported by the CDF collaboration ||, through the ratio Rtb, 



Bit -> bW) 



0.94 ± 0.27(stat) ± 0.13(syst) . 



(1) 



Z q = d , s , b B(t^qW) 



Assuming three generations, this yields 



V tb 



0.97 ±0.15 ±0.07 , 



(2) 
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which is consistent with unity but within experimental errors also consistent with a value 
which is considerably less than unity, namely at 95% C.L. one gets \Vtb\ > 0.58. This 
measurement is expected to improve significantly in future. A precision of 5|Vf&|/|Vf&| ~ 
12% is projected at the Fermilab Tevatron with an integrated luminosity of 2(/6)~ 1 , 
expected to be collected at the turn of this century ||. Eventually, \V t b\ will be measured 
in experiments at the linear collider(s) with a precision 5|14&|/|^4&| = 0(1 — 2)% from the 
anticipated accuracy of 5T(t)/T(t) ~ 1% on the top quark decay width 0. However, it 
will be difficult in the foreseeable future to get quantitatively useful information on Vtd 
and V ts from direct decays of the top quark, both due to the anticipated small branching 
ratios involving these matrix elements, 

B(t -> sW) = 0(1(T 3 ), B(t -> dW) = 0(1(T 4 ), (3) 

and, more importantly, due to the (present) low efficiency of tagging light-quark jets. This 
is somewhat discomforting as the direct determination of the CKM matrix elements in 
top quark decays and their inferred values from FCNC processes, such as the ones from 
B decays being discussed in these lectures, would have given very stringent constraints on 
possible new physics or perhaps would have established its existence. It is likely that the 
FCNC processes in B (and to a lesser extent in K) decays will remain the major source 
of information on Vtd and Vt s . We will discuss the present quantitative determinations 
of these matrix elements and their possible improved measurements at the forthcoming 
B facilities, such as the B factories, HERA-B, and the hadron colliders (Tevatron and 
LHC), using rare B decays, AM^ and AM S . 

The weak interaction phase responsible for CP violation in the SM resides dominantly 
in the matrix elements V t d and V u b- This is manifest in the Wolfenstein parameterization 
H of the CKM matrix. Hence B decays and mixings involving one or both of these matrix 
elements are potentially the most promising means to measure CP violation. Since the 
information on the CP violating phase is rather sparse, essentially confined at present to 
the decay Kl — > ttvt, the CP violating asymmetries in B decays will be very welcome 
and perhaps decisive input in testing the CKM paradigm for CP violation. We shall 
give a profile of these CP violating asymmetries in B decays and the underlying CKM 
unitarity triangle based on fitting the present data and will discuss measurements at 
future facilities which will go a long way in reducing the present uncertainties in the CKM 
parameter space. These experiments (and the anticipated theoretical progress) have the 
possibility of putting the quark flavour physics at a comparable precision level as the 
present electroweak physics in the post-LEPI era. 

This writeup is organized as follows: In section 2, we introduce the CKM matrix and 
the unitarity triangle (s) using the Wolfenstein parametrization for this matrix. In section 
3, we discuss the dominant B decay modes, which determine the bulk quantities such 
as the semileptonic branching ratio Bsl(B), the average number of charmed particles 
per B decay (n c ), and the individual B hadron lifetimes. The present determination of 
the matrix elements \V c b\ and \V u b\ from semileptonic B decays is also discussed in this 
section, using the HQET methods for the former. In section 4, we take up the discussion 
of the electromagnetic penguins and rare B decays in the SM and make comparison with 
data in terms of the branching ratio and the photon energy spectrum. This measurement 
determines the ratio of the CKM matrix elements |V^|/|V^,| which we quantify. The 
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CKM-suppressed radiative rare decays B — > + 7 and the corresponding exclusive 
decays B — > (p, 07)7 are discussed in section 5. Their role in determining the CKM matrix 
element \V t d\ (equivalently the CKM-Wolfenstein parameters p and rf) is reviewed. The 
success of this proposal depends in a crucial way on reliable calculations of the so-called 
long distance (LD) contributions and we discuss some existing estimates of the same. In 
this section, we also take up the FCNC semileptonic decay B — > X s £ + £~ in the SM model, 
discussing first the QCD-improved rates and distributions from the short-distance (SD) 
contribution, including leading power corrections. A quantitative understanding of these 
decays requires reliable estimates of the LD and non-perturbative effects which we also 
discuss. In section 6, we give an update of the CKM matrix and the unitarity triangle 
(UT), taking into account the present measurements and theoretical estimates in a number 
of B decays and |e|, the CP violating parameter in the kaon sector. The constraints on the 
CKM parameters from the present LEP bound on AM, are also analyzed. In section 7 we 
briefly discuss some representative CP-violating asymmetries in B decays and summarize 
their expected ranges and correlations in the SM. We conclude with a brief summary in 



section 8. Some of the topics discussed here have also been reviewed in [1C 



2 CKM Matrix and the Unitarity Triangle 

We start by discussing the flavour changing transitions in the SM. Since QCD is manifestly 
flavour-diagonal, the only possibility of FCNC transitions is in the electroweak sector. 
Writing in terms of the physical boson (W^, Z®, A^) and fermion fields, it is easy to 
show that the neutral current part of the standard electroweak model is also manifestly 
flavour-diagonal. Denoting the quarks and leptons by fi(i = 1...6), the neutral current in 
the SM is given by: 



sin 9 W cos 9w 



[hi ~ Qsm 2 6 w ) . + eA^Qi 



fi, (4) 



where = (1 — Isj/^h) with I 3 = +1/2 for Ui and z/j and —1/2 for di and charged 
leptons £i, and Qi is the electric charge of the fermion fi in units of the electron charge, 
i.e., Q e = +1. The electroweak mixing angle in J^ c , denoted by 9w, has its origin 
in the diagonalization of the gauge boson mass matrix, and it has the usual definition 



cos 9w = 92/ y gl + with the electric charge defined as e = g 2 sin 6w Concerning 
the Higgs Yukawa couplings - a potential source of FCNC transitions in general - it is 
known that the unitary transformations which diagonalize the quark mass matrix also 
diagonalize the Higgs Yukawa sector in the SM. This is most easily seen by writing the 
Yukawa sector of the SM Lagrangian, which after spontaneous symmetry breaking has 
the form 



"Yukawa 



■\//; <',/,' + daMtdjR + luMijijRl ( 1 + — ) /'•<•• • ( 5 ) 



where the absence of the neutrino mass matrix is conspicuous and represents the SM 
choice of treating the neutrino massless (equivalently, the absence of the right-handed 
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neutrinos vm). In the basis where the quark masses are diagonal, this takes the form 

£ Yukawa = ~ E ( l + ~) > ( 6 ) 

where H is the Higgs field and vq is the Higgs vacuum expectation value. This manifest 
flavour diagonal form of £ Yukawa ^ n g enera l is n °t maintained in multi-Higgs models and 
one has to impose discrete symmetries on the Higgs and fermion fields to forbid FCNC 
couplings in ^Yukawa' as emphasized by Glashow and Weinberg quite some time ago 
P~H] . The absence of such couplings in the SM owes itself to the choice of a single Higgs 



doublet. 

The charged current J^ c in the SM, which couples to the W^, is 

/ d 

J T = /o ■ n 5 ' *)i ^CKM 



M y/2svciO w 



s | , (7) 

where V^km = V^ p V^ own ^ is a (3x3) unitary matrix in flavour space, first written down by 
Kobayashi and Maskawa in 1973 ||. The matrices V^ p and y^ own diagonalize the up-type 
and down type quark mass matrices, respectively. The matrix Vqkm is a generalization 
of the Cabibbo rotation for the three-quark-flavour (u, d, s) case, invented to keep 
the universality of weak interactions, which took the form of a (2 x 2) matrix by the 
inclusion of c-quark with the GIM construction Q , and is called the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix. There are no FCNC transitions in the SM at the tree level by 
construction. They are induced by higher order CC transitions and the resulting FCNC 
amplitudes are determined by the masses of the intermediate quarks, i.e. they reflect the 
flavour dependence of the Higgs Yukawa couplings, weighted with the appropriate CKM 
prefactors. 

The charged current in the SM has a (V — A) structure, hence it violates P and C 
maximally, conserves the electric charge and the lepton- and baryon-number separately, 
but otherwise there are no restrictions on it except that ^ckM^CKM = 1- I n general, 
C cc violates CP due to the possibility of a non-trivial phase in Vqkm- 
Symbolically the matrix V^KM can be written as: 



/ V u d V us 



v ckm 




V c d v cs 
V v td v ts 

For quantitative discussions we need a parametrization of the CKM matrix. The original 
parametrization due to Kobayashi and Maskawa || was constructed from the rotation 
matrices in the flavour space involving the angles Qi (i = 1, 2, 3) and a phase S, 

Vkm = #2 3 (#3, 5)Ri2{0i, 0)R 2 3(0 2 , 0), (9) 

where < 0, < ir/2, < 5 < 2tt, and Rij(9,4>) denotes a unitary rotation in the 
plane by the angle 9 and the phase 0. The resulting representation is: 

-S]C 3 -S1S3 
%M = I s i c 2 C1C2C3 - s 2 s 3 e tS C1C2S3 + s 2 c 3 e' t5 \ , (10) 
C1S2C3 + c 2 s 3 e tS C1S2S3 - c 2 c 3 e lS 
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with Cj = cos6i,Si = sin^j. This reduces to the usual Cabibbo form for 6 2 = 63 = 0, 
with the angle #1, identified (up to a sign) with the Cabibbo angle. In the PDG review 
12fl , however, another parametrization is advocated which differs from Vj^m in assigning 



the complex phases (dominantly) to the (1,3) and (3,1) matrix elements of Vqkm- An 
approximate but very useful form of the matrix Vqkm is due to Wolfenstein ||: 





^Wolfenstein = | ~^ 1 - ^X AX \ . (11) 

-p-irj) 

with A = sin 0c ~ 0.221. Like other representations, Vwolfenstein h as a l so three real 
parameters called A, X and p, and a phase rj. Since we shall be making extensive use of 
this parametrization, we write some relations involving the matrix elements of interest in 
this representation: 

iKfcl . rz \V t 



a Vp 2 + v 2 , \rn = ^(i-p) 2 + v 2 , (12) 



v cb \ v " " \v cbl 

7? n/ 1 

1 , (13) 



\V td \ /(l-p) 2 + ^ \v t 



ts 



\v ub \ V p 2 + v 2 \v cb 

and the dominant phases are: 



3(Kb) = %(V td ) = -AX\ (14) 



It should be recalled that the Wolfenstein parameterization given in Eq. (|iT|) is an approx- 
imation and in certain situations in the future it may become mandatory to specify the 
matrix by taking into account the dropped terms in 0(A 4 ) in Vv^ if enste j n . For the present 
experimental and theoretical accuracy, the representation (|ll|) is entirely adequate and 
we shall restrict ourselves to this form. Further discussions on this point and suggestions 
on improved treatment to include higher order terms in A can be seen in [13|]. 



2.1 The CKM unitarity triangles 

The CKM matrix elements obey unitarity constraints, which state that any pair of rows, 
or any pair of columns, of the CKM matrix are orthogonal. This leads to six orthogonality 
conditions which can be depicted as six triangles in the complex plane of the CKM 
parameter space ||14|| . The constraint stemming from the orthogonality condition on the 
first and third row of V^kM' 

V ud V td + V us V t : + V ub V t l = (15) 

has received considerable attention. Since, as discussed in the introduction, present mea- 
surements are consistent with V ud — 1, V tb — 1 and V t * s ~ —V cb , the unitarity relation flT5| ) 
simplifies to: 

V ub + V td = V us V cb , (16) 
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which can be conveniently depicted as a triangle relation in the complex plane, as shown 
in Fig. [l|. We shall refer to it as the unitarity triangle (UT). Thus, knowing the sides 
of the UT, the three angles of this triangle a, (3 and 7 are determined. These angles are 
all related to the Kobayashi-Maskawa phase 5 (equivalently the phase #i 3 in Vprjc or 
the phase rj in Vwolfenstein)> and they can, in principle, be independently measured in 
various CP-violating B decays. Restricting to the Wolfenstein representation in which 
the dominant phases reside in the (13) and (31) matrix elements, these angles are defined 
as follows: 



V ub V u \ 2 V (r ] 2 + p 2 -p) 



sin 2a = arg 
sin 2(3 = arg 



sin 2 7 = arg | ) = 2 . (17) 



v: b v; d ) (p 2 + v 2 )((i-p) 2 + v 2 ) 

_VtA 27/(1 -p) 

, v;J (i-p) 2 + v 2 ' 

Vub \ __ 2r]p 

v:J~ P 2 +v 2 



Some estimates of the angles a, (3 and 7, and hence CP asymmetries in B decays, can 
be obtained at present by constraining the parameters of the CKM matrix p and 77. 
Conversely, knowing the CP asymmetries, the parameters p and rj of the CKM matrix 
can be determined. 

As already stated, the matrix elements V cb and V u b are known from the CC B decays. 
With more data and improved theory (in particular for the so-called heavy-to-light tran- 
sitions b — > uivg) one would be able to determine these matrix elements rather precisely. 
The matrix element Vtd can, in principle, be determined from the rare decays b — > d + 7, 
b — > d + / + /~, b — > d + vv (and the corresponding exclusive decays), and B^-B^ mix- 
ing. The mass difference AM^ already provides a first measurement of \V t d\- This set of 
measurements, which involves decay rates and mixing frequencies but not CP-violating 
asymmetries, then provides another way of determining the triangle, namely by measuring 
its sides. The lengths of these sides in the Wolfenstein approximation are 

Rt = yl^l- (18) 

A Vcb 

The CP asymmetries in B decays related to the angles a, (3 and 7 and the sides of the 
unitarity triangle obey the geometric relations: 

sin (3 sin(a + 7) sin (3 

t>i, 



sin a sin a sin(a + (3) 

sin 7 sin(a + /3) sin 7 

K t = : -r- f — — • l iy J 

sin a sma sm(7 + (3) 

By measuring both the sides and the angles, the UT will be overconstrained which is 
one of the principal goals of the current and forthcoming experiments in B physics. Be- 
fore leaving this topic, we note that including the 0(A 5 ) terms in the imaginary part 
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Figure 1: The unitarity triangle. The angles a, j3 and 7 can be measured via CP violation 
in the B system and the sides from the CC- and FCNC-induced B decays. 

of the CKM matrix, an additional phase emerges in the matrix element combination 
arg(— V^Vcb/V^Vcd)- This phase, being equal to X 2 r] is bounded from the CKM fits to 
be less than 0.025 and hence very small. However, showing that this CP violating phase 
is indeed small is both a test of the SM and belongs on the agenda of the forthcoming 
experiments in B decays. 

Obviously, there exists a large number of different parametrizations of the CKM ma- 
trix. However, since the phases of the quark fields are unphysical quantities, the different 
parametrizations, emerging from specific choices of these phases, must all be equivalent. 
The parametrization independent quantities are the absolute values of the matrix elements 
\Vij\ (hence also the angles of the unitarity triangles) and the area of the unitarity trian- 
gles, which is the same for all six triangles and is an invariant measure of CP violation. 
This can be expressed as 

area[A(CKM)] = ^AX e r] (20) 

for the Wolfenstein parametrizations of Vckm- The Jarlskog invariant denoted by the 
symbol J [15] is twice this area, which in the standard model is typically of O(10~ 5 ). It 



is being debated if the intrinsic smallness of J in the standard model is a serious problem 
in explaining the measured baryon asymmetry of the universe (BAU), whose quantitative 
measure is the ratio of the baryon number density to entropy density, 

A B ^, (21) 



with Ab = (0.6 — 1) x 10~ 10 fll2 |. Relating the Jarlskog constant to AB is a profound 
problem and a topic of intense theoretical research. It is an interesting question if B 
decays will lead to some helpful clues towards understanding this connection. Recent 



studies indicate that baryogenesis in the SM at the electroweak scale is unlikely [16] due 
to the LEP constraints on the Higgs boson mass. Very probably, the baryon number 
violation (as well as the lepton number violation) is governed by the physics at the grand 
unification scale, which has then little direct influence on B decays. Experiments in B 
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physics will, however, provide an answer to the question if additional CP violating phases 
in the flavour changing sector are present. These experiments, together with the searches 
of CP violation in the flavour diagonal sector such as the electric dipole moment of the 
neutron, will determine the effective low energy theory of CP violation. 



3 Dominant B Decays in the Standard Model 

We now turn to the mainstream B physics and discuss the dominant decay rates which 
determine the lifetimes of the B hadrons tb, the semileptonic branching ratios Bsl, and 
the charm quark multiplicity in B decays (n c ) - a quantity which has become an important 
ingredient in understanding the semileptonic branching ratios. 

The effective lowest-order weak interaction Hamiltonian can be expressed in terms of 
J^ c , introduced earlier, 

Hw = ih ( j z Cj ^ cc + H > w 

where Gp is the Fermi coupling constant. The calculational framework that is used is 
QCD and we concentrate first on perturbative QCD improvements of the decay rates and 
distributions in B decays. The leading order (in a s ) perturbative QCD improvements us- 
ing TCw have been worked out in semileptonic processes in |L7] - p3| , which are modeled 



on the electromagnetic radiative corrections in the decay of the /z-lepton |24| . For the non- 
leptonic decays, perturbative QCD corrections are calculated in the effective Hamiltonian 
approach using the renormalization group techniques |25j — [^8| . The underlying theoreti- 



cal framework and its numerous applications in weak decays of the K and B mesons have 



been reviewed in a comprehensive paper by Buchalla, Buras and Lautenbacher p9[ , to 
which we refer for details and confine ourselves here to some selected topics. 

Apart from these perturbative QCD improvements, resulting in the so-called QCD- 
improved quark-parton model, one could also improve the quark-parton model itself by 
including power corrections in 1/mq. The method that is used in discussing such cor- 
rections is based on the heavy quark limit of QCD which allows one to do a systematic 
expansion of decay amplitudes in 1/niQ, where wig 3> Aqcd, and Aqcd is the QCD 
scale parameter which is typically of 0(200 MeV) |12| . This technique f3(| - has 
the satisfying feature that the parton model for heavy quark inclusive decays emerges as 
the leading term in the expansion of the decay amplitudes. These methods can also be 
applied to calculate the energy-momentum spectra of the decay products except in the 
end-point region, where the heavy quark expansion breaks down. Here, one has at present 
little choice other than smearing the (theoretical and experimental) spectra with weight 
functions to make meaningful comparison or modeling the non-perturbative effects. We 
shall return to the discussions of these topics later. 



3.1 Inclusive semileptonic decay rates of the B hadrons 

We start with the assumption that the inclusive decays of B hadrons can be modeled 
on the QCD-improved quark model decays. More specifically, while calculating rates, we 
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shall be equating the partial and total decay rates of the B hadrons to the corresponding 
expressions obtained in the parton model, relying on the heavy quark expansion [[30] - 
3§: 



T(B^X) = T(b^ x) + 0(l/m 2 b ) . (23) 
For b quark semileptonic decays, one has two partonic CC transitions: 

b — ► cfv t , (24) 

There exists a close analogy between the b quark decays and /x decay, fi~ — > e~i? e v> fl , 
with the identification: 



[b,(c,u),vi,t ]<->[fi ,e ,^ e ,^]. (25) 

This analogy holds also at the one loop level; 0(a) QED corrections to \i~ decay and 
0(a s ) QCD corrections to b semileptonic decays are related by simply replacing JT7L ITR [H| 



a 



1 s 4 

-a s Tr ^Xi\i = -a s 

6 i=i 6 



(26) 



where Aj are the Gell-Mann SU(3) matrices, and a s is the lowest order QCD effective 
coupling constant, 

12vr , x 

ots = ^r— , (27) 

33-2^)111(^3=^-) 
QCD 

where nf is the number of effective quarks. The semileptonic decay rates can then be 
read off the expression for the 0(a) radiatively corrected /i-decay rate pi |. The rates for 
b — >■ (u, c)£u£ decays, setting me = m Vl = 0, are given by the expression: 



^sdb — ► (u,c)£v e ) = T f(ri) 



1 - 



2a.JmT 



71 



9{r% 



(2? 



with Tq being the normalization factor in the lowest-order rate 



To 



192tt ; 



ni 



Ti = mi/m b (i = u, b), and 



f(r) 



^ + 8r D 



24r lnr. 



(29) 



(30) 



The function g(r) has the normalization g(0) = tt 2 — ^, and numerically g(0.3) ~ 2.51, 
relevant for the b — > u and b — > c transitions, respectively 0, [18|, pLQfl - With Aqqq — 200 
MeV and n/ = 5, this gives about (15)% corrections to the semileptonic decay widths 
involving £ — e,/i : reducing T S l compared to the lowest order result T^l = T f(r). The 
corresponding decrease in the decay width for the semileptonic decay b —>■ ctv t is obtained 
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by an expression very similar to the above one in which the r-mass effects are included 
in the phase space and in the QCD corrections. 



r(6 -> ctu t ) = T P(x c , x T , 0) 



i + 



0) 



(31) 



where P(xi,X2,xs) is the well known three-body phase space factor given for arbitrary 
masses Xi = rrii/mb in |J7|. The function g(xi,X2,xs) has been calculated for arbitrary 



arguments in |21j in terms of a one-dimensional integral. The functions P(xi,0,0) and 
g(x c , 0, 0) go over to the functions /(r) and (— )g(r), respectively, given above for the 
massless lepton case. The numerical values for g(x c , x T , 0) and g(x c , 0, 0) are tabulated in 
38fl . For the default value x c = 0.3, one has g(x c ,x T ,0) = —2.08, yielding about a 12 % 



decrease in T(b — > ctu t ) compared to T^l(b — > cru T ) as a result of the leading order QCD 



corrections ||21|| . For more modern calculations of the decay rate Tsiib — > cru T ), see |22| . 



3.2 Inclusive non-leptonic decay rates of the B hadrons 

The dominant CC-induced non-leptonic and semileptonic decays of B hadrons are gov- 
erned by the effective Lagrangian, 



C 



eff 



-e,fj,,T 



+ h.c. 



(32) 



and we have just discussed the 0(a s ) renormalization effects to the matrix elements of 
the semileptonic piece in £ e //- Here 0\ and 2 are the colour-octet and colour-singlet 
four- Fermi operators, respectively (a and (3 are colour indices), 



01 = (d a Up) L (cf3b a ) L , 

2 = {d a u a ) L (cpbj3) L , 



(33) 



and qi = 1/2(1 — 75) denotes a left-handed quark field. The operators 0[ are related to 
the corresponding fields Oj by the relacement d —> s. The octet-octet (Oi) and singlet- 
singlet (O2) operators emerge due to a single gluon exchange between the weak current 
lines (quark fields) and follow from the colour charge matrix (T^) algebra: 



T ik T ji - —^-jydikSji + -8u8 jk 



(34) 



Here, N c = 3 for QCD. The Wilson coefficients Cj(/i) are calculated at the scale /x = my/ 
and then scaled down to the scale typical for B decays, /i = 0(mt,), using the renormal- 
ization group equations, which brings to the fore the influence of strong interactions on 
the dynamics of weak non-leptonic decays. Without QCD corrections, the two Wilson co- 
efficients have the values C\(mw) = 0, C^irnvi) = 1- Since the operators 0\ and O2 mix 
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under QCD renormalization, it is convenient to introduce the operators 0± = (02±C*i)/2 
having the Wilson coefficients C± which renormalize multiplicatively [25]. The results are 

a s (m w ) - a s {fi) ^ ] f ^ /V 



now known to two-loop accuracy [28 
C±{n) = L±(ii) 



An 



2/?o V^ ) fa 



+ ots{m w ) 



47T 



(35) 



where the multiplicative factor in this expression represents the solution of the RG equa- 
tions in the leading order QCD p5|, 



a s {M ] 



w , 



ld± 



a s {n) 



(36) 



and the exponents have the values d + = 7+ /(2/3 ), cL = 7^ /(2/3 ). The quantities 7± 
are the coefficients of the anomalous dimensions involving the operators 0± (and 0' ± ), 



with 



(o) 
7+ 



4, 7 



(0) 



7?> 



+ 7« 



-14 



8 
9 



(37) 



(3? 



in the naive dimensional regularization (NDR) scheme, i.e., with anticommuting 75. The 
f3i are the first two coefficients of the QCD /5-function, and they have the values 



/5 Q = 11 - /3i = 102 - yn ; . 



(39) 



Finally, the functions _B± are the matching conditions obtained by demanding the equality 
of the matrix elements of the effective Lagrangian calculated at the scale fi = mw and in 
the full theory (i.e., SM) up to terms of 0(a s (m^/)). They have the values: 



Ba 



±B 



2N r . 



(40) 



The constant B and the two-loop anomalous dimension 7^ are both regularization- 
scheme dependent. In the NDR scheme one has B = 11. Following [p8 |, we define a 
scheme-independent quantity i?±, 



#4 



B, 



(0) , (1) 

7± T± 



2(3 V 7 f A) 



(41) 



in terms of which the Wilson coefficients read 

oc s (m w ) - «s(^) 



C±(/x) = L±((i) 



1 + 



An 



47T 



(42) 



In this form all the scheme-dependence resides in the coefficients B± which is to be can- 
celled by the scheme-dependence of the matrix elements of the corresponding operators. 
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In addition to the decays b — > c + ud, b — > c + us and b — > c + £ise, which are described 



by the effective Lagrangian fl32|) , there are other decays involving the CC transitions 
b — > uX, b — > (c, w) + cs and b — > (c, w) +cd, which are not included in this Lagrangian. In 
a systematic treatment involving QCD renormalization, one has to enlarge the operator 
basis to include these transitions and the so-called penguin operators. We shall return to 
a discussion of this part of the Lagrangian later in these lectures as we discuss rare B- 
decays, where the operator basis will be enlarged and the corresponding Wilson coefficients 
calculated in the leading logarithmic approximation. 

We now discuss the semileptonic branching ratio *BgL for the B mesons and to be 
specific will consider the case t = e, //. This branching ratio is to a large extent free of 
the CKM matrix element uncertainties but requires a QCD-improved calculation of the 
inclusive decay rates, Tg^ , discussed above, and 



SL 



F(B -> Xcv e 

r to t(£) 



(43) 



with 



r tot(£) = E T(B^X£v e )+T(B->X c X) + T(B^X c - c X) 

+ T(B -> X U X) + T(B) (Penguins) . (44) 

In the spirit of the parton model, we shall equate T(B — > X C X) = T(b — > cud) + T(b —>■ 
cus), noting that the so-called ^-annihilation and H^-exchange two-body decays are 
expected to be small in inclusive B decays. This will be quantified later as we discuss the 
lifetime differences among B hadrons which arise from the matrix elements of the operators 
representing these contributions. The corrections for the decay widths T(b — > cud) and 
T(b — > cus) are identical neglecting m u and m s , and so their contributions can be described 
by similar functions. The resulting next-to-leading order QCD corrected sum can be 
expressed as: 



r(6 — > cud) + T(b — > cus) 



x 



+ 



r P(x c ,o,o) 

2L(u)l + L(ut + ^ Mw) - asM ^,,2 , rr.,2 



2tt 



2a s (fi) /3 



(2L(u) 2 + R + + L(u)_R- 
3 



~ L(u)_) 2 c n (x c ) + -(L(u) + + L(u)„) 2 c 22 (x c ) 



3tt 

+ - L(u) 2 _)(c 12 (x c , u) - 12 In—' 

= 3T ri(u)J(x c ,u) , 

with r](a) representing the leading order QCD corrections. The scheme independent R± 
come from the NLO renormalization group evolution and are given by 



(45) 



10863 - 127871/ + 80^ 
6(33 - 2n f ) 2 
15021 - 1530n f + 80n 2 f 



7 



3(33 - 2n 



f) 



(46) 
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For n f = 5, R + = 6473/3174, R_ = -9371/1587. Note that the leading dependence of 
L{n)± on the renormalization scale /i is canceled to 0(a s ) by the explicit /i-dependence 
in the a; s -correction terms. Virtual gluon and Bremsstrahlung corrections to the matrix 
elements of four fermion operators are contained in the mass dependent functions c 



X 



The analytic expressions for the functions Cn(x), c^ix), 022(2;) are given in where also 
their numerical values are tabulated. Lumping together all the perturbative and finite 
charm quark corrections in a multiplicative factor A c (m b , x c , a s (mz)), the perturbatively 
corrected decay width can be expressed as: 



r(6 — > cud) + T(b — > cus) 

For the central values of the parameters used here (m b = 4.8 GeV, x 
Oi s {rnz) = 0.117), the QCD corrections lead to an enhancement [35 

0.17 



3r P(x c , 0, 0) [1 + A c (x c , m b , a s {m z ))} . 

0.3, /i 



A c (m bl x c , a s {m z )) 
Out of this, the bulk is contributed by the leading log factor 



r?0) - 1 



2Li + L_ 



1 = 0.10 . 



(47) 
m b and 

(48) 
(49) 



Next, we equate T(B — > X c5 ) = T(b — > ccs) + T(b — * ccd) and discuss the perturbative 
QCD corrections to the decay width T(b — > ccs) and T(b — > ccd). Neglecting and m s , 
an assumption which has been found to be valid to a high accuracy in the corrections 
in the two decay widths are identical and the result can be written in close analogy with 
the ones for the decay widths T(b — > cus) discussed above. 



r(6 — > ccs) + T(b — > ccd) = T P(x c ,x c ,x s ) 



x 



+ 



a s (M w ) - a s (/i) 



2a s (/i) 



3tt 
1 



2vr 

3 3 

-(L(/i) + - L(fi)_) 2 k u (x c , fi) + + L_) 2 k 22 (x c 



- -(L 2 + -L 2 _)(k l2 (x c ) 



12 In 



/i , 



(50) 



The functions kij(m b , x c , a s (mz)) have been calculated and their numerical values are 
tabulated in ||40|| . Again, lumping together all the perturbative and finite charm quark 
corrections in a multiplicative factor A cc (mb, x c , a s (mz j), the perturbatively corrected 
decay width can be expressed as: 



T(b -> ccs) = 3T P(x c , x c , x s ) [1 + A cc (x c , m b , a s {m z ))\ . 



(51) 



With the values of the parameters used above, the QCD corrections lead to the following 



enhancement 39 



A cc (m 6 , x c , a s (m z )) = 0.37. (52) 

This is by far the largest correction to the inclusive rates we have discussed so far. Using 
pole quark masses and the renormalization scale /i = m b , one gets 

r(6 -> ccs)(NLO) 



91: 



T(b^ccs)(LO) 



1.32 ±0.07 



(53) 
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The NLO corrections go in the right direction in bringing theoretical estimates closer to 
the experimental value for the semileptonic branching ratio. However, this will also lead 
to enhanced charmed quark multiplicity (n c ) in B decays, as discussed a little later. 

The CKM-suppressed and penguin transitions contribute at a smaller rate to T^ Q ^(B). 
They are of two kinds: 

• T(B — > X u + X), which is suppressed due to the CKM matrix element \V u b\, with 
the rate depending on |V^| 2 , and 

• r (B) (Penguin): The so-called penguin transitions b — > s + X, where X = cc and 
X = g (QCD penguins), X = 7 (electromagnetic penguins), X = t+t'^vv (elec- 
troweak penguins). 

There are also transitions involving b — > d + X, as well as a host of other rare decays, 
which can all be neglected. The dominant contributions in the SM add up to [ 1Q| : 

F(B ^X u + X) + r (B) (Penguins) ~ 1.25 x KT 2 r , (54) 

and hence not of much consequence for the semileptonic branching ratio or the B hadron 
lifetime estimates. 



3.3 Power corrections in Ysl{B) and Vnl(B) 

Before we discuss the numerical results for £>sl, we include the 0(l/mf) power corrections 
in the inclusive partonic decay widths. They constitute the first non-trivial corrections to 
the parton model results and have been calculated using the operator product expansion 
techniques |3(|]- [[35)1 . 

In HQET, the 6-quark field is represented by a four-velocity-dependent field, denoted 
by b v (x). To first order in l/m&, the 6-quark field in QCD b(x) and the HQET- field b v (x) 
are related through: 



b{x) 



-im^v.x 



2rrih 



bJx) 



The QCD Lagrangian for the b quark in HQET in this order is: 
^HQET 



b v iv. lpb v + b v — b v — Zbb v — y~ b v + (J 



2mi 



4m b 



mi 



(55) 



(56) 



where Zb is a renormalization factor, with Z\,{\i = mj) = 1 and p = -D M 7 M , with 
being the covariant derivative. The operator b v (i p) 2 b v /2mb is not renormalized due to 
the symmetries of HQET. (In technical jargon, this is termed as a consequence of the 
reparametrization invariance of £ H Q ET .) With this Lagrangian, it has been shown in 
| 30| - that in the heavy quark expansion in order (1/m 2 ), the hadronic corrections 
can be expressed in terms of two matrix elements 



(B^\b v {iD)\\B^) = 2m BW Ai, 
(B^\b v 9 -a^F» v b v \BM) = 2d BW m BW A 2 , 



(57) 
(58) 
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where F^ u is the gluonic field strength tensor, and the constants d B (*) have the value 3 and 
— 1 for B and B*, respectively. The constant A2 can be related to the hyperfine splitting 
in the B mesons, which gives: 



A, 



\m\, -m 2 B ) = 0.12 GeV 2 . 



(59) 



The other quantity Ai is the average kinetic energy of the b quark inside a B meson and 
has been estimated in various ways, using the QCD sum rule approach [IT], the virial 
theorem [g|], lattice QCD |43| and data jPJ. A range 

Ai = -(0.5 ± 0.3) GeV 2 (60) 

is compatible with most estimates. (For a recent compilation of Ai estimates, see [|I3"| .) 
Taking into account these corrections, the semileptonic and non-leptonic decay rates of a 
B meson B — > Xiug and B — > X C X can be written as [ |3~T| , |3"2] ] : 



r(B — Xc&//) = r(°V(r c ) 



1 



2« S K) 



7T 



-g{r c 



1 + 



Ai 
2m 2 b 



+ 



3A 2 
2m 2 b 



6(1 



r c ) 4 A 2 



f(r c 



a. 



1 1 s 



m 



b / 

(61) 



and 



r(s — > x c x) 



3r (o) 



r/(/i)J(/i) 1 + 



Ai 
2mg 



+ 



3A 2 
2mg 



6(1 



r c ) 4 A 2 



/(r c 



m 



- L + (/i) 2 -L_(/i) ; 



4(1 -r c ) 3 A 2 2 a 

— k + Cla 



/(r c 



mf mi 



(62) 



where the product r/(yu)J(/i) denotes the NLO corrected result for the partonic decay 
discussed above in (|45|), to which Eq. (^) reduces in the limit Ai = A 2 = 0. 

The decay rates depend on the quark masses, which unlike lepton masses, do not 
appear as poles in the 5*-matrix nor do the quarks exist as asymptotic states. They 
are parameters of an interacting theory and hence subject to renormalization effects. 
Consequently, they require a regularization scheme, such as the MS scheme, and a scale, 
where they are normalized, to become well-defined quantities. For example, the quark 
masses in the so-called MS scheme and the pole masses (OS scheme) are related in the 



leading order 45 



ojjno)_ 

(37T) + - 



(63) 



In HQET, quark masses can be expressed in terms of the heavy meson masses itlm 
and the parameters Aj, A 2 and a quantity called A, where 



This yields 



m b - m c 



- Ai + g?mA 2 

m M = m Q + A h 

2m Q 

Ai + 3A 2/ 1 



•m B - m D + 



1 



+ 0(- 



(64) 



(65) 
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and the quark mass differences can then be calculated knowing Ai and A2, giving (mj — 
m c ) = (3.4 ± 0.03 ± 0.03) GeV 4(|. This difference, which determines the inclusive 
rates and shape of the lepton energy spectrum in semileptonic decays, has also been 
determined from an analysis of the experimental lepton energy spectrum in B decays, 
yielding (m b — m c ) = 3.39 ± 0.01 GeV for the pole masses [44], in excellent agreement 
with the QCD sum rule based estimates. 



3.4 Numerical estimates of Bsl(B) and 



\ n c) 

The theoretical framework described in the previous section can now be used to predict 
two important quantities in B decays Bsl{B) and (n c ), which have been measured. Con- 
cerning Bsl{B), there is some discrepancy between the two set of experiments performed 
at the T(AS) and at the Z° resonance, although it must be stressed that these experiments 
measure a different mixture of B hadrons. The present measurements give: 

B SL (B) = (10.37 ± ±0.30)% at T (AS) [47], 

B SL {B) = (11.11 ±0.23)% atZ @, 

(n c ) = 1.16 ±0.05 at T(4S) H 

(n c ) = 1.23 ±0.07 at Z° (66) 

where the number for (n c ) at the Z° is from the ALEPH collaboration. We use the follow- 
ing average in which the error on Bsl{B) is inflated to bridge the gap in the experimental 



measurements 50 



B SL (B) = (10.90 ± ±0.46)% 
(n c ) = 1.18 ±0.04. 



(67) 



The theoretical predictions for these quantities have been updated by Bagan et al. |3 



and more recently by Neubert and Sachrajda IpO] , using the same theoretical input. We 
shall use here the numerical results from |50| where the following ranges of parameters 
have been used: 



m b (pole) = 4.8 ± 0.2 GeV; a s {m z ) = 0.117 ± 0.004, m b /2 < ji < 2m b 



(6* 



and 0.25 < m c /m b < 0.33. Here m b is the pole mass defined to one-loop order in 
perturbation theory. At order 1 /ml in the heavy quark expansion, non-perturbative 
effects are described by the parameter A2, as the dependence on the parameter Ai cancels 
out in calculating Bsl(B) and (n c ). This analysis leads to the following values using the 
pole masses (OS scheme) |5( 



B S L 

(n c ) 



(12.0 ±1.0)% (for// 
1.20 =f 0.06 (for /j, = 



■ m b ), (10.9 ±1.0)% (for fi = m b /2) 
m b ), 1.21 =f 0.06 (for p = m b /2). 



(69) 



One could also use, following Bagan et al. |$9|], the MS scheme and the results in this 
scheme are as follows f|50H 



B SL (MS) = (10.9 ±0.9)% (for (i = m b ), (10.3 ± 0.9)% (for /i = m b /2) 
(n c ) = 1.25 =F 0.05 {for fi = m b ), 1.24 =f 0.06 (for ft = m b /2), 



(70) 
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The numbers in the MS scheme correspond to using the two-loop anomalous dimension 
matrix in the running of the quark masses and the errors from various sources have been 
added in quadrature. The estimates ( |69D and (|70|) show that both (n c ) and Bsl are 
scheme-dependent; in addition Bsl also depends on the scale /i. A comparison of the 
theoretical estimates in the OS scheme (eqs. fl69|)) and data on (n c ) and Bsl is shown 
in Fig. [| Given the parametric dependence on the scale /i and the ratio m c /mb, the 
agreement between theory and experiment is reasonably good. In the MS scheme, the 
semileptonic branching ratio is generally smaller and (n c ) somewhat higher (the two are 
anti-correlated). To make more precise predictions, one has to calculate the missing 0(af) 

1.4 



1.3 
c° 1.2 



1.1 



1 

8 9 10 11 12 13 14 

B SL (%) 

Figure 2: The charm content n c vs. Bsl in B decays, shown as a function of the quark 
mass ratio m c /mt, and the renormalization scale n/rrib in the OS scheme. The data show 
the experimental averages from the T(4S) (LE) and Z° (HE) measurements. (Figure 
taken from pD|). 

corrections, which, as the experience has it, will considerably reduce the scale dependence. 
However, equally important is to reduce the present theoretical uncertainty in the ratio 
m c /mfc. Here, precise data on the lepton and hadron energy distribution in semileptonic 
B decays will help. So, while there is certainly much room for improvement, it is fair to 
conclude that within existing uncertainties the current theoretical estimates for Bsl and 
(n c ) in the SM do not disagree significantly from the corresponding experimental values. 

3.5 B-Hadron Lifetimes in the Standard Model 

A matter closely related to the semileptonic branching ratios is that of the individual 
B hadron lifetimes. The QCD-improved spectator model gives almost equal lifetimes. 
Power corrections will split the 5-baryon lifetime from those of B%, B^ and B s . However, 
first estimates of these differences are at the few per cent level PH| . The experimental 
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situation has been summarized as of summer 1996 in [|5"T|]: 

' {B ' -1.04 ±0.04; = 0.98 ± 0.05; ^4 = 0.78 ± 0.04 . (71) 



r(B d ) * ' r(B d ) ' ' r(B d 

The subject of exclusive B hadron lifetimes has received renewed theoretical attention 
lately 0, 0, in which the possibly enhanced roles of the four-Fermion operators 
involving baryonic states (as compared to the mesonic state) has been studied. We recall 



that such operators enter at 0(1 /ml) in the heavy quark expansion discussed above pi 



In this order, there are four such operators, which using the notation of |)(J, can be 
expressed as: 

O q v _ A = (b Llll q L )(q L Yb L ), 

O q s _ P = (b R q L )(q L b R ), 

V_ A = (b LllJl t a q L )(q Ll H a b L ), 

Ts-p = (b R t a q L )(q L t a b R ), (72) 

where t a are generators of colour SU(3). The matrix elements of these operators between 
various B-meson and A^-baryons are in general different and this contribution will thus 
split the decay widths of the various B hadrons. In general, the operators ( [72] ) introduce 
eight new parameters corresponding to the matrix elements of these operators. In the 



large- N c limit, however, it has been argued in [50] that the B- mesonic matrix elements of 
the operators (B q \O v _ A \B q ) and (B q \O q s _ P \B q ) are the dominant ones. While accurate 
numerical estimates require a precise knowledge of these matrix elements, one expects that 
they give rise typically to the spectator-type effects (using the parton model language): 



r spec _ , 27t/b 2 
r tot m B 



(73) 



with fs of order 200 MeV. In the case of A b baryons, one can use the heavy quark spin 
symmetry to derive two relations involving the operators given above taken between the 
A b states. The problem is then reduced to the estimate of two matrix elements which in 
are taken to be the following: 

1 (A b \O v _ A \A b ) = -^r&), (74) 



2m A / U| 48 K B, 

and 

(A 6 |a^|A 6 > = -B(A b \O q v _ A \A b ) , (75) 

The operator Oy-A is a linear combination of the operators Ty_ A and Oy_ A introduced 
earlier, Oy_ A = 2T v _ A + 30y_ A , following from colour matrix algebra |50|| , and r(A b /B, 



is the ratio of the squares of the wave functions which can be expressed in terms of the 
probability of finding a light quark at the location of a b quark inside A& baryon and the 
B meson, i.e. 



V \< b \ 2 



(76) 
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One expects B = 1 in the valence-quark approximation. However, the ratio r(A b /B q ) has 
a large uncertainty on it, ranging from r(A b /B q ) ~ 0.5 in the non-relativistic quark model 
|53fl to r(A b /B q ) = 1.8 ± 0.5 if one uses the ratio of the spin splittings between E 6 and E£ 



baryons and B and B* mesons, as advocated by Rosner pi and using the preliminary 



data from DELPHI, m(S*) - m(E b ) = (56 ± 16) MeV f5|. 

Using the ball-park estimates that B and r(A b /B q ) are both of order unity yields for the 



lifetime ratio T(A b )/r(B c i) > 0.9 150], significantly larger than the present world average. 
Reliable estimates of these constants can be got, in principle, using lattice-QCD and QCD 
sum rules. Very recently, QCD sum rules have been used to estimate (A b \O v _ A \A b ) and 
B, yielding (A b \6 q v _ A \A b ) = (0.4 - 1.2) x 10~ 3 GeV 3 and B = 1.0 |56|. This corresponds 
to the parameter r(A b /B q ) having a value in the range r(A b /B q ) ~ 0.1 — 0.3, much too 
small to explain the observed lifetime difference. We mention here the possibility of linear 
power corrections in the inclusive decay rates, which are not encountered in the explicit 
power corrections discussed above but may enter via the breakdown of the parton-hadron 



duality. Phenomenological parametrizations presented in [[57] in support of such a scenario 
are interesting but not persuasive. One must conclude that the lifetime ratio r{A b )/T{B c i) 
remains a puzzle. New and improved measurements are needed, which we trust will be 
forthcoming from HERA-B and the Tevatron experiments in not-too-distant a future. 

Before leaving this section, we mention that from a theoretical point of view one 
expects measurable lifetime differences between the two mass eigenstates of the B° s -B® 
complex |3T|] . Recently, leading order perturbative (0(AQQjy/m b )) and power corrections 
(0(m s /m b )) to the differences in the decay rates Tf and T| of the two mass eigenstates have 
been analyzed in [ 5S| . The perturbative corrections go along very much the same lines as 



discussed earlier. Power corrections bring in the four-quark operators already mentioned. 
Quantitative estimate of Ar requires the knowledge of non-perturbative quantities, bag 
factors called B and Bs, involving the expectation values of the operators (Oy_ A ) and 
(O s s _ P ), and the pseudoscalar meson coupling constant f Bs - The resulting expression for 
the ratio (AT/T) Bs can be expressed as : 

(^) R . |aB + 6Bs + cl (_i|_) 2 , (77) 

where the constants B and Bs are the mentioned bag factors and the coefficients a, b and c 
depend on the parameters such as m b and /x; c incorporates the explicit l/m b corrections. 
For the choice B = Bs = 1 (corresponding to the vacuum insertion approximation), 
m b = 4.8 GeV, /i = m b and f Bs = 210 MeV, one gets a = 0.009, b = 0.211, c = -0.065, 



yielding [58 



(-) s> = 0.155. (78) 

This difference is large enough to be measured in the forthcoming experiments. If accu- 
rately measured, AT Bs has the potential of providing an alternative estimate of the mass 
difference in the B° s - B® complex, AM# s , as the ratio ATbJ AMb s , as opposed to the 
mass difference itself, does not depend on fs s ■ However, there is still some dependence in 
this ratio on the unknown bag constants and the coefficients. The present determination 



of this ratio is: AT B J AM Bs = (5.6 ± 2.6) x 10 [pq| , which is in need of substantial 
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improvement if it has to make any impact on the CKM phenomenology. We shall return 
to the estimates of AM S and a discussion of the related issues later. 



3.6 Determination of \V c b\ and \V u b\ 

The CKM matrix element V c b can be obtained from semileptonic decays of B mesons. We 
shall restrict ourselves to the methods based on HQET |59| , |60[1 to calculate the exclusive 
semileptonic decay rates and use the heavy quark expansion to estimate the inclusive rates. 
Concerning exclusive decays, we recall that in the heavy quark limit (m& — > oo), it has 
been observed that all hadronic form factors in the semileptonic decays B — ► (D, D*)iv^ 
can be expressed in terms of a single function, the Isgur-Wise function It has been 
shown that the HQET-based method works best for B — > D*lv decays, since these are 
unaffected by 1/mg corrections f. 
B -> D*iv t is 



Using HQET, the differential decay rate in 



dT(B -> D*k 
du 



48vT 



(m B - m D *) 2 vr? D *r] 2 A \J J 2 - l(u + if 



(79) 



x [1 + 



4a; l-2ur + r 2 
u + 1 (1-r) 2 J 



where r = mp* /rriB, uj = v ■ v' (v and t> ' are the four- velocities of the B and -D* meson, 
respectively), and i]a is the short- distance correction to the axial vector form factor. In 
the leading logarithmic approximation, this was calculated by Shifman and Voloshin some 
time ago - the so-called hybrid anomalous dimension pi] . In the absence of any power 
corrections, £(u = 1) = 1. Estimating the size of the 0{l/m\) and 0(l/m 2 c ) corrections 
to the Isgur-Wise function, £ (u) has generated some lively theoretical debate |^SJ, EG] but 



it seems that a convergence has now emerged on their magnitude. We take [46]: 



£(1) = 1 + S(l/m 2 ) = 0.945 ± 0.025 . 



10) 



The quantity and its counterpart for the vector current matrix element renormaliza- 
tion, f]y, have now been calculated in the complete next-to-leading order by Czarnecki 
getting 



r] A = 0.960 ±0.007, 
T] V = 1.022 ±0.004. 



The error on .F(l) is now dominated by the power corrections in £(1), yielding |66 

=£. VA = 0.907 ±0.026 . 



(81) 



(82) 



Since the rate is zero at the kinematic point u = 1, one uses data for uo > 1 and an 
extrapolation procedure to determine £(1)|V^|. As the range of accessible energies in the 
decay B — > D*£u is rather small (1 < uo < 1.5), the extrapolation to the symmetry point 
can be done using a Taylor expansion around u = 1, 



= ^(1) 1 - p 2 {uo - 1) + c{uo - l) 2 + ... 



*3) 



20 



The present experimental input from the exclusive semileptonic channels is based on 
the data from CLEO, ARGUS, ALEPH, DELPHI and OPAL, which is summarized by 
Gibbons at the Warsaw conference [67|. He obtains 



cb\ 



0.0357 ±0.0020 ±0.0014 . 



34) 



Using from Eq. (p2[) , gives the following value: 

\V cb \ = 0.0393 ± 0.0021 (expt) ± 0.0015(curv) ± 0.0017 (th), 



35) 



where the error from the curvature of the Isgur-Wise function has also been indicated. 
Combining the errors quadratically gives 



w, 



cb\ 



0.0393 ± 0.0028 



(86) 



A value of [V^l has also been obtained from the inclusive semileptonic B decays using 
heavy quark expansion. The inclusive analysis has the advantage of having very small 
statistical error. However, as discussed previously, there is about 2a discrepancy between 
the semileptonic branching ratios at the T(4S) and in Z° decays. Using an averaged value 
for the semileptonic decay width from these two sets of measurements and inflating the 
error as before to take into account the disagreement leads to a value ||47|| : 



\V cb \ = 0.0398 ± 0.0008 (expt) ± 0.004 (th) 



*7) 



where the theoretical error estimate (±10%) has been taken from Neubert f4"6fl . For further 
discussion of these matters we refer to |65], £46|. The agreement in the values of \V c b\ 
obtained from the exclusive and inclusive semileptonic decays is remarkably good. This 
can be taken as a quantitative check of the notion of quark-hadron duality in semileptonic 
decays. We shall use the following values for \V c b\ and the Wolfenstein parameter A in 
the CKM fits discussed later: 



w, 



eft 



0.0393 ± 0.0028 



A = 0.81 ±0.058 



The knowledge of the CKM matrix element ratio iKib/Kbl is based on the analysis 
of the end-point lepton energy spectrum in semileptonic decays B — > X u lvi and the 
measurement of the exclusive semileptonic decays B — > (tt,p)£i>i reported by the CLEO 
collaboration [47]. As noted in |68 |, the inclusive measurements suffer from a large ex- 
trapolation factor from the measured end-point rate to the total branching ratio, which is 
model dependent. The exclusive measurements allow a discrimination among a number of 
models all of which were previously allowed from the inclusive decay analysis alone. 
It is difficult to combine the exclusive and inclusive measurements to get a combined de- 
termination of [Kifel/I V cb \. However, it has been noted that the disfavoured models in the 
context of the exclusive decays are also those which introduce a larger theoretical disper- 
sion in the interpretation of the inclusive B — > X u ivg data. Excluding them from further 
consideration, measurements in both the inclusive and exclusive modes are compatible 



with 57 



ub 



0.08 ±0.016 



39) 
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This gives 

yp 2 + r] 2 = 0.363 ± 0.073 . (90) 

We summarize this section by observing that the bulk properties of B decays are largely 
accounted for in the standard model. On the theoretical front, parton model estimates 
of the earlier epoch have been replaced by theoretically better founded calculations with 
controlled errors. In particular, methods based on HQET and heavy quark expansion 
have led to a quantitative determination of |V c f,| at ±7% accuracy, which makes it after 
\V u d\ and \V US \, the third best measured CKM matrix element. The matrix element \V u t\ 
has still large uncertainties (±20%) and there is every need to reduce this, as this error is 
one of the two handicaps at present in testing the unitarity of the CKM matrix precisely. 
The quantities Bsl-, (n c ), and the individual S-hadron lifetimes are now in reasonable 
agreement with data. A completely quantitative comparison requires the missing NLL 
corrections and in the case of lifetime differences better evaluations of the matrix elements 
of four-quark operators, which we hope will be forthcoming. Finally, we stress that it will 
be very helpful to measure the semileptonic branching ratios Bsl for the baryons. 
With the lifetimes of the B hadrons now well measured, such a measurement would allow 
to compare T S L(B d ), r 5i (_B ± ) and r Si (A b ), to check the pattern of power corrections in 
semileptonic decays. 

4 Electromagnetic Penguins and Rare B Decays in 
the Standard Model 

We now discuss the FCNC transitions SM which in the SM are induced through the 
exchange of W bosons in loop diagrams. We shall discuss representative examples from 
several such transitions involving B decays, starting with the decay B — > X s + 7, which 
has been measured by CLEO ||69|| . This was preceded by the measurement of the exclusive 



decay B -+ K* + 7 |7g|: 

B(B^X s + -f) = (2.32 ± 0.57 ± 0.35) x 10" 4 , (91) 
B{B^K* + j) = (4.2 ±0.8 ±0.6) x 10" 5 , (92) 

yielding an exclusive-to-inclusive ratio: 

T(B -> K* + 7) 

R -'- r { B^x, + ',) =m±6 - 8)% - (93) 

These decay rates test the SM and the models for decay form factors and we shall study 
them quantitatively. 

The leading contribution to b — > s + 7 arises at one-loop from the so-called penguin 
diagrams and the matrix element in the lowest order can be written as: 

M(b^s + 7) = -| JL £ V lb V*F 2 (xi) q^sa^(m b R + m s L)b , (94) 

where Xi = m 2 /m 2 v , and e M are, respectively, the photon four-momentum and po- 
larization vector, the sum is over the quarks, u, c, and t, and V%j are the CKM matrix 



22 



elements. The (modified) Inami-Lim function i^Xj) derived from the (1-loop) penguin 
diagrams is given by fTD: 



24(x 



6x(3a;-2)loga;-(x-l)(8a; 2 + 5a;-7) , (95) 



where in writing the expression for F 2 (xi) above we have left out a constant from the 
function derived by Inami and Lim, since on using the unitarity constraint these sum to 
zero. It is instructive to write the unitarity constraint for the decays B — > X s + 7 in full: 

V tb V* + V cb V; s + V ub V* s = . (96) 

Now, since the last term in this sum is completely negligible compared to the others (by 
direct experimental measurements), one could set it to zero enabling us to express the 
one-loop electromagnetic penguin amplitude as follows: 

M(b -> s + 7) = -| ^2 A * - F ^ x ^ Q^ v sa^(m b R + m s L)b . (97) 

The GIM mechanism || is manifest in this amplitude and the CKM-matrix element 
dependence is factorized in A* = Vt b V t *. The measurement of the branching ratio for 
B — > X s + 7 can then be readily interpreted in terms of the CKM-matrix element product 
At/|V c fe| or equivalently |Vt s |/|Vd,|. In the approximation we are using (i.e., setting X u = 
0), this is equivalent to measuring \V CS \. For a quantitative determination of |Vt s |/|K*|, 
however, QCD radiative corrections have to be computed and the contribution of the 
so-called long-distance effects estimated. We proceed to discuss them below. 



4.1 The effective Hamiltonian for B — > X s j 

The appropriate framework to incorporate QCD corrections is that of an effective theory 
obtained by integrating out the heavy degrees of freedom, which in the present context 
are the top quark and bosons. This effective theory is an expansion in 1/m^ and 
involves a tower of increasing higher dimensional operators built from the quark fields 
(u, d, s, c, b), photon, gluons and leptons. The presence of the top quark and of the W ± 
bosons is reflected through the effective coefficients of these operators which become 
functions of their masses. The operator basis depends on the underlying theory and in 
these lectures we shall concentrate on the standard model. The basis that we shall use is 
restricted to dimension-6 operators and the operators which vanish on using the equations 
of motion are not included. The effective Hamiltonian T~C e ff given below covers not only 
the decay b — ► s + 7, in which we are principally interested in this section, but also other 
processes such as b — > s + g and b — > s + qq. 

It is to be expected in general that due to QCD corrections, which induce operator- 
mixing, additional contributions with different CKM pre-factors have to be included in 
the amplitudes. Thus, QCD effects alter the CKM-matrix element dependence of the 
decay rates for both B — > X s + 7 and (more importantly) B — > Xd + 7. However, with the 
help of the unitarity condition given above, the CKM matrix dependence in the effective 
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Hamiltonian incorporating the QCD corrections for the decays B — > X s + 7 factorizes, 
and one can write this Hamiltonian as 

AC 8 

W (6 _> s + 7 ) = -^^^Q^^C/x), (98) 

V2 4=1 

where the operator basis is chosen to be (here \i and v are Lorentz indices and a and j3 
are colour indices) 

01 = {sLcTfpb^fapfcLp), (99) 

02 = (sLa^b Lf3 )(c L ^C La ), (100) 

03 = feaT/Ac*) E fe/^M, (101) 

04 = (*£a7„M E (fe/^Va), (102) 

05 = feaT/Aa) E (^7^), (103) 

06 = {sLalM E (S^y^Jto), (104) 

°7 = T^misLaa^bi^F^, (105) 

°7 = T^m^s^^fex^F^, (106) 



16tt 2 

08 = ^m 6 (s iQ T^^&^)G a ^, (107) 

° 8 = J§^ m ^R» T >^)G a>1 ^ (108) 

where e and g s are the electromagnetic and the strong coupling constants, and F^ v and 
G^j, denote the electromagnetic and the gluonic field strength tensors, respectively. We 
call attention to the explicit mass factors in 7 {0' 7 ) and 8 (0g), which will undergo 
renormalization just as the Wilson coefficients. The dominant contributions in the radia- 
tive decays B — > X s + 7 arise from the operators 02, 07 and 0$, whereas the operators 
03, 06 get coefficients through operator mixing only, which numerically are negligible. 



Historically, the anomalous dimension matrix was calculated in a truncated basis [72 



and this basis is still often used for the sake of ease in calculating the real and virtual 
corrections, though as we discuss below, now the complete anomalous dimension matrix 
is available [73| . 

The perturbative QCD corrections to the decay rate T(B — > X s + 7) have two distinct 
contributions: 

• Corrections to the Wilson coefficients Cj(/i), calculated with the help of the renor- 
malization group equation, whose solution requires the knowledge of the anomalous 
dimension matrix in a given order in a s . 



^Note that in addition to the penguins with the it-quark intermediate state there are also non- 
factorizing contributions due to the operators (wl q 7 m 6l q )(sl^7 ai ui / ^), which like the it-quark contribution 
to the 1-loop electromagnetic penguins are proportional to the CKM- factor X u = V us V* b , and hence are 
consistently set to zero. 
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• Corrections to the matrix elements of the operators Oi entering through the effective 
Hamiltonian at the scale /i = 0(m&). 

The anomalous dimension matrix is needed in order to use the renormalization group and 
sum up large logarithms, i.e., terms like a™(mw) log m (mb/M), where M = m t or rrtw 
and m < n (with n = 0, 1, 2, ...). Until recently, only the leading logarithmic corrections 
(m = n) have been calculated systematically in the complete basis given above [[73|] . Very 
recently, the next-to-leading order anomalous dimension has also been calculated and 



reported this summer by Misiak at the Warsaw conference j74 . 

Next-to-leading order corrections to the matrix elements are now available completely. 
They are of two kinds: 

• QCD Bremsstrahlung corrections b — > 57 + g, which are needed both to cancel the 
infrared divergences in the decay rate for B — > X s + 7 and in obtaining a non-trivial 
QCD contribution to the photon energy spectrum in the inclusive decay B — > X s + r y. 

• Next-to- leading order virtual corrections to the matrix elements in the decay b — > 
s + 7. 



The Bremsstrahlung corrections were calculated in [75[ - [|77| in the truncated basis and 
last year also in the complete operator basis ||78|1 , which have been checked in [79|. The 
higher order matching conditions, i.e., Ci(mw), are known up to the desired accuracy, 
i.e., up to 0(a s (Mw) terms [|80|. The next-to- leading order virtual corrections have been 
completed by Greub, Hurth and Wyler recently []8T|. We discuss the presently available 
pieces in the SM calculation of BiB —>■ X s + 7) in the NLO accuracy. 

We recall that the Wilson coefficients obey the renormalization group equation 



9 n , \ 9 



a 



,9) = iji(g)Cj 



fj 



(109) 



A* / \ /i 

The QCD beta function /3(g) has been defined earlier and 7(g) is the anomalous dimension 
matrix, which, to leading logarithmic accuracy, is given by 

,2 



7(0) = 7o 



g- 



16tt 2 



(110) 



Here 70 is a 8 x 8 matrix given in [7 
given at the scale M w and read |7T] 



The non-zero initial conditions in the SM are 



c 2 


= 1 


C 7 (M W ) 


1 

= — X 
2 


C 8 (M W ) 


1 

= — X 
2 



2x 2 /3 + 53712 - 7/12 3x 2 /2 - x 

{x - l) 4 
3x/2 



\nx 



x 



x 2 /4-5x/4- 1/2 
(x- l) 3 



+ 



In a; 



[x 



(111) 
(112) 

(113) 



and x = m^/Mvy- Also, for subsequent discussion it is useful to define two effective Wilson 
coefficients C7 (p) and Cg (fi) [f^]: 



Cf = 



c 7 - 



3 



c fi 



(114) 
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The numerical values for the Wilson coefficients at \x = My/ ("Matching Conditions") 
and at three other scales \x = 10.0 GeV, 5.0 GeV and 10.0 GeV can, for example, be seen 



in JTOf and will not be given here. In LO, one gets [32|: 

B{B -> X s + 7) = (2.8 ± 0.8) x 10~ 4 , 

reflecting the parametric uncertainties of the underlying framework. 

Now, we discuss the explicit 0(a s ) improvement of the decay rate. The real and 
virtual 0(a s ) corrections to the matrix element for b — > s + 7 at the scale \x ~ m;, form 
a well-defined gauge invariant albeit scheme- dependent set of corrections. This scheme 
dependence will be cancelled against the one in the anomalous dimension 7' 1 ', as discussed 
in the context of the dominant contributions to the non-leptonic decays of the B hadron 
earlier. The results presented here correspond to the NDR scheme. 

The Bremsstrahlung corrections in b —>■ 57 + g, calculated in [|75j - [77[| and ]75 |, were 
aimed at getting a non-trivial photon energy spectrum at the partonic level. In these 
papers, the virtual corrections to b — > 57 in 0(a s ) were included only partially by taking 
into account those virtual diagrams which are needed to cancel the infrared singularities 
(and also the collinear ones in the limit m s — ► 0) generated by the Bremsstrahlung 
diagram. The emphasis was on deriving the photon energy spectrum in B — > X s + 7 
away from the end-point x 7 — > 1 and the Sudakov-improved photon energy spectrum in 
the region x 1 — > 1. The left-out virtual diagrams, however, do contribute to the overall 
decay rate in B — > X s + 7. Recently, these virtual correction have been evaluated in |8l| , 
neglecting the small contributions from the operators O^-Oq. The additional contribution 
reduces substantially the scale dependence of the leading order (or partial next-to-leading 
order) decay width T(B —>■ X s + 7), which previously was found to be substantial and 
constituted a good fraction of the theoretical uncertainty in the inclusive decay rate 



82, 84, 78 



Concentrating on the dominant operators O2, Oj and Og, 
next-to-leading order correction to the matrix element part in b - 
as follows: 

M = M 2 + M 7 + M 8 



the contribution of the 
-> s + 7 can be expressed 

(115) 



and the various terms (including appropriate counter term contributions) can be summa- 



rized as f81| : 



Mk 



(Sj\0 7 \b) tr 



with 




1 mb A 
2 log h r 2 



(116) 



;ii7) 



Kr 2 



243 
+ 

+ 

+ 



{-833 + 1447rV /2 

1728 - 180tt 2 - 1296C(3) + (1296 - 324tt 2 )L + 108L 2 + 36L E 
648 + 72?r 2 + (432 - 216tt 2 )L + 36L E 
-54 - 84vr 2 + 1092L - 756L : 
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Sr 2 



16tt ( r 



45 - 3vr 2 + 9L + 9L 2 ] z + [-3vr 2 + 9L 2 ] z 2 + [28 - 12L] z 3 }ll9) 



Here, 3?r2 and Qr2 denote the real and the imaginary part of r2, respectively, z = (m c /mb) 2 
and L = log(z). 

The real and virtual corrections associated with the operator 7 , calculated in J7^, 
76| , 78| can be combined into a modified matrix element for 6 — > 57, in such a way that 
its square reproduces the result derived in these papers. This modified matrix element 
My od reads [0: 



M n 7 



mod 



with 



(sj\O r \b) tree \ l + ^U logy + r 7 



r 7 = - (4 - Ti 
7 9 V 



Finally, the result for J\A$ is [81 



a. 



m b 



M 8 = {si\0 7 \b) tree -a- 4 log -± + r 8 



Arc 



with 



-33 + 2tt 2 - 6ztt 



(120) 
(121) 

(122) 
(123) 



32 4 
* = - T , r 8 = -- 

With the results given above, one can write down the amplitude M.(b — > S7) by 
summing the various contributions already mentioned. Since the relevant scale for a b 
quark decay is expected to be \i ~ m b , the matrix elements of the operators may be 
expanded around \i = m b up to order 0(a s ) and the next-to- leading order result can be 
written as: 

4G F A t 



M(b — > S7) 



v/2 



L>( S7 |0 7 (m b )|6)i re 



with 



D = cv e ^) + ^^(a (0)e// ' 



/i 



47T 



*(0)e//„ 



(124) 



(125) 



where the quantities 7$*^ — i% + 8 ^7 are the entries of the (effective) leading order 
anomalous dimension matrix and the quantities ii and r, are given for i = 2, 7, 8 in eqs. 
CTJTT8T ), ( |T2l ) and Q), respectively. The first term, C 7 e# (/i), on the r.h.s. of Eq. ( |T25| ) 
has to be taken up to next-to-leading logarithmic precision in order to get the full next- 
to- leading logarithmic result, whereas it is sufficient to use the leading logarithmic values 
of the other Wilson coefficients in Eq. ( |1 25| ) . 

The decay width Y mrt which follows from Ai(b — > sj) in Eq. ( p.24j ) reads 



•vuirt _ m b,pole G 2 F \ 2 a em 

~ 32^ 



F\D\ 



(126) 



where the terms of 0(a 2 ) in \D\ 2 have been discarded. The factor F in Eq. ( |126| ) is 



' m b (/j, = m b ) 



a s {m b ) 



7T 



(127) 
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and its origin lies in the explicit presence of mj in the operator O7. To get the inclusive 
decay width for b — > s r y(g), also the Bremsstrahlung corrections (except the part already 
absorbed above) must be added. The contribution of the operators 2 and O-j was cal- 



culated already in [|75|]. As pointed out by Buras et al. p2| , the explicit logarithms of 
the form a s (mj,) log(m b //i) in Eq. ( |125| ) are cancelled by the /i-dependence of C^ 6 ^ (n). 
Therefore, the scale dependence is significantly reduced by including the virtual correc- 
tions completely to this order. This is shown in Fig. |3] (solid curves), which yields an 
error of ±6% on BiB — > X s + 7) varying /1 in the range m^jl < fi < Im^. We recall that 
in the LO calculations, this scale dependence was ±20%. The CLEO result is shown as 
dashed lines. The two other curves (dashed-dotted) represent more stringent assumption 
on the /^-dependence and we refer to [HI]] for further details. 
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Figure 3: Branching ratio for b — > s / y(g) calculated in |8T[] with the parameters |T4 S |/|K 
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4.8 GeV and m c /m{, = 0.29. The different curves are explained in 



4.2 Estimating long-distance effects in B — > X s + 7 

In order to get the complete amplitude for B — > X s + 7 one has to include also the effects 
of the long-distance contributions, which arise from the matrix elements of the four-quark 
operators in 7i e //, (X s j\Oi\B). It is usual to invoke the hypothesis of factorization, which 
is then combined with the additional assumption of vector meson dominance, involving 
the decays B — > + X s — >7 + X s , where Vi = J/ip,?P', ■■■ H^] - |89| . One has to 



ensure that the resulting amplitude A4(B — > X s + 7) remains manifestly gauge invariant. 
In practice, this amounts to discarding the longitudinal polarization contribution in the 
non-leptonic decays B — > (J/ip, ip', ...)+X s , which in fact dominates the decay widths |90|| . 
and keeping only the smaller contribution from the transverse polarization of J/ip,i/)', .... 
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Following WA B3], one can write the decay amplitude as: 



M(b -> sJ/^t = ^agV^ ^^^ a + 7fe)fc&4fa) , (128) 
where is defined as (-0(g) jcy^cIO) = —igii,(q 2 )tt(q)- F° r the decays under consideration 



one needs the value of the BSW coefficient a 2 ||91|| , which has been determined to be 
\a 2 \ = 0.24 ± 0.04 p0| . One also needs to evaluate the coupling constant gv(q 2 ) at the 
point q 2 = 0. From leptonic decays of vector mesons, one gets, however, gy{q 2 = My). 
As noted in [|87|, [88|], using gv(q 2 = 0) = gv(q 2 = My) would substantially overestimate 



the long-distance contribution due to the expected dynamical suppression of the effective 
coupling gv(q 2 ), as one extrapolates to the point q 2 = 0. In fact, such a suppression is 
supported by HERA data on photoproduction of J/ip. Including all the (cc) resonances 
and the short distance contribution Msd, the two-body decay amplitude {b — > sj) can 
be written as 

M(b^ S7 ) = - e ^V tb V: s [^- 2 m b D^) - a 2 | £ 4^]^(1 + 7b) 6^ , (129) 

where ipi represents the following vector cc resonant states: i/;(lS), i[>(2S), ^(3770), 
^(4040), ^(4160), and ^(4415), and D is the function given earlier. Taking this estimate 
as giving the right order of magnitude for the long-distance contribution, Deshpande et 
al. |[7| conclude that such long-distance effects can be as large as 10%. Other estimates, 
in particular by Golowich and Pakvasa p8[ , lead to an even smaller long-distance con- 



tribution. Clearly, one can not argue very conclusively if such estimates are completely 
quantitative due to the assumptions involved. In future, one could improve them by using 
data from HERA on elastic J/ip—, and ^'-photoproduction to get gj/^(0) and 5y(0) di- 
rectly, reducing at least the extrapolation uncertainties involved in the presently adopted 
procedure of extracting these coupling constants from the leptonic decay widths of each 
state and extrapolating to the point q 2 = using an Ansatz. In conclusion, LD effects in 
B — > X s + 7 are dynamically suppressed. 

4.3 Estimates of B(B X s + 7) in the Standard Model 

In the quantitative estimates of the SM branching ratio B(B — > X s + 7) given below 
we have neglected the LD-contributions. It is theoretically preferable to calculate this 
quantity in terms of the semileptonic decay branching ratio 

B{B - X sl ) = [ r(jB ~7 + Xs) r B(B - Xivi) , (130) 

1 SL 

where, the leading-order QCD corrected Tsl has been given earlier. The leading order 
power corrections in the heavy quark expansion, discussed in the context of the semilep- 
tonic decay rate, are identical in the inclusive decay rates for B — > X s + 7 and B — > X£i>£, 
entering in the numerator and denominator in the square bracket, respectively, and hence 



drop out |3Q, 31 



The error on the branching ratio B[B — ► A s 7) comes from the following sources. 



29 



1. Am t and A/x: The present value of m t is m t = 175 ± 9 GeV ||, which is usually 
interpreted as the pole mass. With this the running top quark mass in the MS 
scheme is rn t = 166 ± 9 GeV. This leads to an error of ±4% in B(B — ► X s j). The 
combined error on Am t and Afj, is about ±9% as can be seen in Fig. || 

2. Errors from the extrinsic parameters (A(m&), A(a s (rriz)), and A(BRsl), the exper- 
imental uncertainty on the semileptonic branching ratio): This gives an uncertainty 
of ±12% on B[B — > X s ±7) as estimated in |7^| , of which half is due to the assumed 
uncertainty A(a s {m z )) = 0.006. 

As mentioned already, the Wilson coefficient has been calculated in the next-to- 
leading order |74|. The NLO corrections are found to be small but the numerical difference 



is still being worked out. Replacing Cy by its leading log value yields the branching ratio 

B(B -f X s + 7) = (3.20 ± 0.30 ± 0.38) x 10~ 4 (131) 

where the first error comes from the combined error on Am t and A/i, as can be seen in 
Fig. H and the second from the extrinsic source. Combining the theoretical errors in 
quadrature gives 

B(B -> X s + 7) = (3.20 ± 0.48) x 10" 4 . (132) 

Using the same input, a branching ratio BiB — > X s + 7) = (3.25 ± 0.50) x 10~ 4 has been 
calculated in [85fl. The error on the NLO branching ratio is now ±15%, reduced by a 



factor 2 from the corresponding LO value. The SM branching ratio B(B — > X s ± 7) is 
compatible with the present measurement BiB — > X s + 7) = (2.32 ± 0.67) x 10~ 4 |^ . 



On its face value, the electroweak penguin rate in the SM is nominally larger than the 
present experimental value, but due to the large errors this difference is not significant. 
Nevertheless, this comparison suggests that there is little room for an additive beyond- 
the-SM contribution, which, for example, is the case in multi-Higgs doublet models. For 
constraints on such models, see Expressed in terms of the CKM matrix element 

ratio, one gets 

1^1 



w, 



0.85±0.12(expt) ±0.08(th), (133) 



cb\ 



which is within errors consistent with unity, as expected from the unitarity of the CKM 
matrix. 

Finally, we note that the ratio Rk* has been calculated in a large number of models. 
Not surprisingly, taken together they give rise to a large dispersion for this quantity. 
However, one should stress that QCD sum rules and models based on quark-hadron 
duality give values which are in good agreement with the CLEO measurements. Some 
representative results are: 



Rk* 


= 0.20 ±0.06 


Rk* 


= 0.17 ±0.05 


Rk* 


= 0.16 ±0.05 


Rk* 


= 0.16 ±0.05 


Rk* 


= 0.13 ±0.03 



(134) 
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4.4 Photon energy spectrum in B — > X s + 7 

The two-body partonic process 6 — > 37 yields a photon energy spectrum which is just a 
discrete line, l/(r)<fT(& — > 57) = 5(1 — x), where the scaled photon energy x is defined as 
Ej = (m 2 — m 2 ) / (2 mj) x. The physical photon energy spectrum is built up by convoluting 
the non-perturbative effects due to the hadronic states involved in the decay and the per- 
turbative QCD corrections, such as b — > s'j+g, which give a characteristic Bremsstrahlung 
spectrum in x in the interval [0, 1] peaking near the end-points, E y — > E™ ax (or x — > 1) 
and Ej — > (or x — > 0), arising from the soft-gluon and soft-photon configurations, re- 
spectively. Near the end-points, one has to improve the spectrum obtained in fixed order 
perturbation theory. This is usually done in the region x — > 1 by isolating and exponen- 
tiating the leading behaviour in a em a s (fj,) m log 2n (l — x) with m < n, where \i is a typical 
momentum in the decay B — > X s + 7. The running of a s is a non-leading effect, but 
as it is characteristic of QCD it modifies the Sudakov-improved end-point photon energy 



spectrum |95|, [96[ compared to its analogue in QED |p7|1 . As long as the s-quark mass is 
non-zero, there is no collinear singularity in the spectrum. However, parts of the spectrum 
have large logarithms of the form a s \og(m 2 /m 2 ), which are important near the end-point 
x — > but their influence persists also in the intermediate photon energy region and they 



have to be resummed [75, |98j. Implementation of non-perturbative effects is at present 
model dependent. 

We shall confine ourselves to the discussion of the photon energy spectrum calculated 
in |75|, [77], |7^| , incorporating the perturbatively computed spectrum, discussed in the 
previous section, and non-perturbative smearing for which a model is used. In this model 
fl9|l , which admittedly is simplistic but has received some theoretical support in the HQET 
approach subsequently |[2| [L00|| , the b quark in B hadron is assumed to have a Gaussian 
distributed Fermi motion determined by a non-perturbative parameter, pp, 

4 — p 2 

<t>{P) = 7= n exp( — -) , p=\p\ (135) 

yffPF PF 2 

with the wave function normalization / °° dpp 2 (f)(p) = 1. The photon energy spectrum 
from the decay of the 5-meson at rest is then given by 



dE 1 



dpp 2 ^p)^W lVl E,) , (136) 



where p ma x is the maximally allowed value of p and is the photon energy spectrum 
from the decay of the 6-quark in flight, having a momentum-dependent mass W(p). This is 
calculated in perturbation theory taking into account the appropriate Sudakov behaviour 
in the E y end-point region at the partonic level. 



An analysis of the CLEO photon energy spectrum has been undertaken in ]78| to 
determine the non-perturbative parameters of this model, namely m^pole) and pp. The 
latter is related to the kinetic energy parameter Ai defined earlier in the HQET approach. 
The experimental errors are still large and the fits result in relatively small x 2 values; 
the minimum, \ 2 is obtained for pp = 450 MeV and mb(pole) = 4.77 GeV. While the 
value of the kinetic energy parameter Ai is at present in a state of flux and hence 
no quantitative conclusions can be drawn, the value of the 6-quark pole mass determined 
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from the photon energy spectrum is in good agreement with theoretical estimates of the 
same, namely m^ipole) = 4.8 ±0.15 GeV [ |101| , |102j| . In Fig. |we have plotted the photon 



energy spectrum normalized to unit area in the interval between 1.95 GeV and 2.95 GeV 
for the parameters which correspond to the minimum x 2 (solid curve) and for another set 
of parameters that lies near the x 2 -boundary defined by x 2 — Xmin + 1- (dashed curve). 



Data from CLEO |69| are also shown. Further details of this analysis can be seen in 

It is a very desirable goal to calculate the photon energy spectrum in B — > X s + 7 in a 
theoretically more robust framework. In this context we note that attempts to calculate 
the photon and lepton energy spectra in the heavy quark expansion method lead to formal 

1Q0| . Near the end-point, the 



expressions which near the end-point are divergent |32|, [ 
energy released for the light quark system in the final state is not of 0{m^) but of the 
order of the parameter A = Mb — ra^ ~ O(Aqqj)). Thus, the expansion parameter is 
no longer l/m 2 , but rather l/m^A = 0(1) and the operator product expansion breaks 
down. This divergent series in the effective theory has to be cleverly resummed and the 
distributions averaged over momentum bins [ |103| |. The resummation allows us to define, 
in principle, an effective non-perturbative shape function [ [L03| , |9~6| , which though can not 



be calculated in the effective theory but one could use this concept advantageously to 
relate the energy spectra in the semileptonic decays B — > X u tvs> and B — > X s + 7. Since 
the perturbative corrections are process-dependent, they have to be included accordingly. 
This programme remains to be implemented in a phenomenological analysis of data. 



p F = 450 MeV, m q = 

p r = 310 MeV, m = 300 MeV 




2.4 2.6 
E ? [GeV] 



Fig. 4e 



Figure 4: Comparison of the normalized photon energy distribution using the CLEO data 
|69| corrected for detector effects and theoretical distributions from JF8 , both normalized 
to unit area in the photon energy interval between 1.95 GeV and 2.95 GeV. The solid 



curve corresponds to the values with the minimum x , ( m q-, 
dashed curve to the values (m q ,p F ) = (300 MeV, 310 MeV). 



p F ) = (0, 450 MeV), and the 



5 Inclusive radiative decays B —> + 7 

The theoretical interest in studying the (CKM-suppressed) inclusive radiative decays B — > 
Xd + 7 lies in the first place in the possibility of determining the parameters of the CKM 
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matrix. With that goal in view, one of the interesting quantities in the decays B — > X d + ^y 
is the end-point photon energy spectrum, which has to be measured requiring that the 
hadronic system X d recoiling against the photon does not contain strange hadrons to 
separate the large- E 1 photons from the decay B — > X s + 7. Assuming that this is feasible, 
one can determine from the ratio of the decay rates B(B —> X d + r y)/B(B — > X s + 7) 
the CKM-Wolfenstein parameters p and 77. This measurement was first proposed in |76| , 
where the photon energy spectra were also worked out. 

In close analogy with the B — > X s + 7 case discussed earlier, the complete set of 
dimension-6 operators relevant for the processes b — > dry and b — > d-yg can be written as: 

AC 8 

n eff (b - d) = Y, c Mdifa), (is?) 

where £j = VjbV* d with j = t,c,u. The operators Oj, j = 1,2, have implicit in them 
CKM factors. In the Wolfenstein parametrization |J, one can express these factors as : 

Z u = A\ 3 (j>-i7i), Zc = -A\\ & = (138) 

We note that all three CKM-angle-dependent quantities £ 3 - are of the same order of mag- 
nitude, 0(A 3 ). It is calculationally convenient to define the operators 0\ and O2 entering 



in Tl e ff(b — > d) as follows [76 



01 = —f{c L prf'b Lo ){d Lo r( li c L p) - —{u L p^b La ){d L(X ^ IJi u L p), 

it « 

2 = —^(c La YbLa){dL^CL/3) ~ («L a 7 M &La) ( rf L/37^/3) > ( 139 ) 

it it 

with the rest of the operators Oj defined like their counterparts Oj in H e ff(b — > s), with 
the obvious replacement s d. With this choice, the matching conditions Cj(mw) and 
the solutions of the RG equations yielding Cj(fi) become identical for the two operator 
bases Oj and Oj. The essential difference between T(B — > X s + 7) and r(5 — > X d + 7) 
lies in the matrix elements of the first two operators 0\ and O2 (in 7i e ff(b s)) and 0\ 
and O2 (in 7i e ff(b —yd)). The branching ratio £>(5 — > X d + 7) in the SM can be written 
as: 

B{B ~^X d + 1 ) = D 1 \ 2 {(1 - pf + rj 2 ~{\- p)D 2 - V D 3 + D 4 }, (140) 

where the functions Di depend on the parameters m t , nib, m c , p, as well as the others we 
discussed in the context of B(B — >• X s + 7). These functions were first calculated in |HJ 
in the leading logarithmic approximation. Recently, these estimates have been improved 
in [104] , making use of the NLO calculations in JS]] discussed in the context of the decay 



B — ► X s + 7 earlier. To get the inclusive branching ratio, the CKM parameters p and 
7] have to be constrained from the unitarity fits. Present data and theory restrict the 
parameters p and rj to lie in the following range (at 95% C.L.) M: 

0.20 < n< 0.52, 
-0.35 < p < 0.35 , (141) 
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which, on using the current lower bound from LEP on the B® - mass difference 
AM S > 9.2 (ps) ( - 1} |7|, restricts p to lie in the range —0.25 < p < 0.35, with r] not 
changed significantly. This is based on assuming £ s = 1.1, where £ s is the S77(3)-breaking 



parameter £ s = f Bs 



Bs d - The preferred CKM-fit values are [|] 
(p, 77) = (0.05, 0.36), 



(142) 



for which one gets 104 



B(B -f X d + 1 ) = 1.62 x 10~ 5 , (143) 

whereas B[B — > X d + 7) = 8.0 x 10~ 6 and 2.8 x 10~ 5 for the other two extremes 
p = 0.35, 77 = 0.50 and p = —i] = —0.25, respectively. In conclusion, we note that 
the functional dependence of B(B — > X d + 7) on the Wolfenstein parameters (p, rf) is 
mathematically different than that of AM S . However, since the no n- factorizing terms 
represented by the coefficients D2 - D± in the expression for B(B — > X d + 7) are nu- 
merically small [104], the resulting constraints from this decay mode and AM d / AM S are 
qualitatively very similar. From the experimental point of view, the situation p < is 
favourable for both these measurements as in this case one expects (relatively) smaller 
values for AM S and larger values for the branching ratio BiB — > X d + 7), as compared 
to the p > case which would yield larger AM S and smaller B(B — > X d + 7). 



5.1 B(B — > V + 7) and constraints on the CKM parameters 

Exclusive radiative B decays B — > V + 7, with V = K*,p,u, are also potentially very 
interesting from the point of view of determining the CKM parameters [p3]| . The extrac- 
tion of these parameters would, however, involve a trustworthy estimate of the SD- and 
LD-contributions in the decay amplitudes. 

The SD-contributions in the exclusive decays (B ± , B°) -> (K**, K*°) +7, (B ± , B°) -> 
(p ± ,p°)+7, B° — > 07 + 7 and the corresponding B s decays, B s — > + 7, and B s — > K*° + j, 
involve the magnetic moment operator Oj and the related one obtained by the obvious 
change s — *> d, O7. The transition form factors governing these decays can be generically 
defined as: 

(V, \\^a^b\B) = ie, vpa ei x V B p° v F*-* v (0). (144) 

Here V is a vector meson with the polarization vector e^ x \ V = p,cu,K* or 0; B is a 
generic 5-meson B^, B° or B s , and ip stands for the field of a light u,d or s quark. The 
vectors pb, Pv and q = Pb — Pv correspond to the 4-momenta of the initial 5-meson and 
the outgoing vector meson and photon, respectively. Keeping only the SD-contribution 
leads to obvious relations among the exclusive decay rates, exemplified here by the decay 



rates for {B ± ,B°) 


-(p ± , 


p°) + 7 and 


(B ± , B )- 
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r((5 ± ,5°) 
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v t 



ts\ 



(145) 



where is a phase-space factor which in all cases is close to 1 and 

Ki = [F s l ^ p / Fg l ^ K ] 2 . The transition form factors Fs are model dependent. Estimates 
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of Fg l ^ K in the QCD sum rule approach are in good agreement with the CLEO data, as 
already discussed. The ratios of the form factors, i.e. Ki, should therefore also be reliably 
calculable in this approach as they depend essentially only on the SU(3)-breaking effects. 

If the SD-amplitudes were the only contributions, the measurements of the CKM- 
suppressed radiative decays (B ± ,B°) — > (p ± ,p°) + 7, B° — > uj + 7 and B s — > K* + 7 
could be used in conjunction with the decays {B^, B°) — > (K , K*°) + / y to determine the 
CKM parameters. The present experimental upper limits on the CKM ratio |V^|/|Vt s | 
from radiative B decays are indeed based on this assumption, yielding at 90% C.L. ||105|| : 



ts 



< 0.45-0.56 



depending on the models used for the SU(3) breaking effects in the form factors [p3 , 94 



The possibility of significant LD-contributions in radiative B decays from the light 
quark intermediate states has been raised in a number of papers [|S6 |. Their amplitudes 
necessarily involve other CKM matrix elements and hence the simple factorization of the 
decay rates in terms of the CKM factors involving |V^| and |V^ S | no longer holds thereby 
invalidating the relation ( [L45| ) given above. As we already discussed, the LD-contributions 
are small in the exclusive decays B — > K* + 7 and so this issue hinges sensitively upon 
the LD-contributions in the CKM-suppressed decays, B^ — > p ±r y and B° — ► (p°,co>)7. 

The LD-contributions in B — > V + 7, induced by the matrix elements of the four- 
Fermion operators 0\ and O2 (likewise 0\ and O2), have been investigated in [|106| , |107|1 
using a technique ||108|| which treats the photon emission from the light quarks in a 
theoretically consistent and model-independent way. This has been combined with the 
light-cone QCD sum rule approach to calculate both the SD and LD — parity conserving 
and parity violating — amplitudes in the decays (B ± , B°) — > (p ± , p/uj) + 7. To illustrate 
this, we concentrate on the B ± decays B ± — > p ± + 7, and take up the neutral B decays 
B° — > p(uS) + 7 at the end. 

The LD-amplitude of the four-Fermion operators 0±, O2 is dominated by the contri- 
bution of the weak annihilation of valence quarks in the B meson and it is color-allowed 
for the decays of charged B ± mesons. Using factorization, the LD-amplitude in the decay 
B — >• p^ + 7 can be written in terms of the form factors F± and F, 



L 

2 ! 



A, 



eG, 



long 



y/2 - au - ua V< • N c 

x{-i[g^(q-p)- I f q v ]-2F 1 L (q 2 
Estimates from the light-cone QCD sum rules give ||107|| : 



+ e 



p aq( i-2F 2 L (q 2 )}. (147) 



F£/F s = (1.25 ± 0.10) x 10" 2 , F^/F s = (1.55 ± 0.10) x 10" 



(148) 



where the errors correspond only to the variation of the Borel parameter in the QCD sum 
rules. Including other possible uncertainties, one expects an accuracy of order 20% for 
the ratios in (|148|) . The parity-conserving and parity- violating amplitudes turn out to be 
numerically close to each other in the QCD sum rule approach, F^ ~ F 2 L = F^, hence 



the ratio of the LD- and the SD- contributions reduces to a number [107] 

y VubV* d 



Along/ A, 



short 



R 



L/S 



v tb v; d 



(149) 
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Using C 2 = 1.10, Ci = -0.235, = -0.306 (corresponding to the scale // = 5 GeV) 
gives: 

pB± ^ ±7 _ 47T 2 m p (C 2 + d/JV.) pB ± ->p ± 'y 

R l/s = = -°- 30 ± °- 07 > ( 15 °) 

which is not small. To get a ball-park estimate of the ratio ^4iong/-4.short, we take the 
central value from the CKM fits, |V^&|/|Kd| — 0.33 0, yielding 

\A long /A short \ B± ^ = |<V P±7 I^717T - 10% • ( 151 ) 

1 \V td Vtb\ 

Thus, the CKM factors suppress the LD- contributions in — ► p ± 7- 

The analogous LD-contributions in the neutral B decays B° —>■ fry and B° — > co>7 are 
expected to be much smaller, a point that has also been noted in the context of the VMD 
and quark model based estimates [|S6j] . The corresponding form factors for the decays 
B° — > p°(c<j)7 are obtained from the ones for the decay B^ — > p ± 7 discussed above by the 
replacement of the light quark charges e u — > e^, which gives the factor —1/2; in addition, 
and more importantly, the LD-contribution to the neutral B decays is colour-suppressed, 
which reflects itself through the replacement of the BSW-coefficient a\ by a 2 ||91|| . This 
yields for the ratio 

r>B°->iry 

4^ = ^^-0.13 ±0.05, (152) 



where the numbers are based on using a 2 /ai = 0.27 ± 0.10 fl90" |. Thus, in this approach 
R-L/s* n — R-l/s*^ 1 = 0-05, which in turn gives 



.4 



long 



A 



short 



< 0.02. (153) 



The above estimate, as well as the one in eq. ( |151|) , should be taken only as indicative 
in view of the approximations made in |106| , |107|1 . That the LD-effects remain small in 
B° — > p7 has been supported in a recent analysis based on the soft-scattering of the on- 
shell hadronic decay products B° — ► p°p° — > fry [ |109|| , though this paper estimates them 
somewhat higher (between 4 — 8%). 

The relations, which follow from the SD-contribution and isospin invariance 

T(B ± -> p ± 1 ) = 2 T(B° -> p° 7 ) = 2 T^ -> oo-y) , (154) 

on including the LD-contributions get modified to 



r(£?± _> p ± 7 ) _ r(i? ± -> P * 7 ) 

2r(5° -> p 7 ) 21(5° -> cj 7 ) 



1 + i? 



1 + 2 • + w ^^fe • (155) 



where -Rl/s = R^/s^ p 1 ■ The ratio 



r(5 ± -> p ± 7 )/2r( J B° -> p 7 )(= r(S ± -> p±7)/2r(5° -> cj 7 )), estimated to lie in the 



36 



range 0.7 - 1.2 in [107|, constrains the Wolfenstein parameters (p, 77), with the dependence 



on p more marked than on 77. 

The ratio of the CKM-suppressed and CKM-allowed decay rates for charged B mesons 
also gets modified due to the LD contributions. Following p6[ , we ignore the LD- 
contributions in T(B — > K*^). The ratio of the decay rates in question can therefore 
be written as: 



T(B ± -> p* 7 ) 

r(5± 




n u \ 2 [{\ - pf + ^ 



2 • R L/B VJ$^* + (RMj^^} , (156) 



Using the central value from the estimates k u = 0.59 ± 0.08 p3 |, we show the ratio ( |1 56|) 



in Fig. |5] as a function of p for rj = 0.2, 0.3, and 0.4. It is seen that the dependence 
of this ratio is weak on 77 but it depends on p rather sensitively. The effect of the LD- 
contributions is modest but not negligible, introducing an uncertainty comparable to the 
0(15%) uncertainty in the overall normalization due to the S'[/(3)-breaking effects in the 
quantity k u . 




Figure 5: Ratio of the CKM-suppressed and CKM-allowed radiative 5-decay rates 
T(B ± -> p ± 7)/r(73 -> K*i) (with B = B ± or B°) as a function of the Wolfenstein 
parameter p, a) with 77 = 0.2 (short-dashed curve), rj = 0.3 (solid curve), and rj = 0.4 
(long-dashed curve). (Figure taken from ||107|| .) 



Neutral S-meson radiative decays are less-prone to the LD-effects as argued above, 
and hence one expects that to a good approximation (say, of O(10%)) the ratio of the 
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decay rates for neutral B meson obtained in the approximation of SD-dominance remains 
valid [0: 

T(B° -y p°7,w7) 2 2 2 

T(B -+ K*>y) = KdX 1(1 ~ p) +V] > (157) 

where this relation holds for each of the two decay modes separately. 

Finally, combining the estimates for the LD- and SD-form factors in ||107|| and |93| , 



respectively, and restricting the Wolfenstein parameters in the range —0.25 < p < 0.35 
and 0.2 < rj < 0.4 from the CKM-fits ||, one gets the following ranges for the branching 
ratios: 

6(5^^7) = (1.5 ±1.1) x 10~ 6 , 
B(B° -> pj) ~ B(B° -> uj) = (0.65 ±0.35) x 10~ 6 , (158) 



where we have used the experimental value for the branching ratio B(B — > K* + 7) |7(J 



adding the errors in quadrature. The large error reflects the poor knowledge of the CKM 
matrix elements and hence experimental determination of these branching ratios will put 
rather stringent constraints on the Wolfenstein parameter p. 

In addition to studying the radiative penguin decays of the B ± and B° mesons dis- 
cussed above, hadron machines such as HERA-B will be in a position to study the corre- 
sponding decays of the B® meson and A& baryon, such as B® — > (ft + 7 and A& — > A + 7, 
which have not been measured so far. Their estimates can be seen in 111 10 . 



5.2 Inclusive decays B -> X s £ + £~ in the SM 

The decays B — > X s £ + £~, with £ = e,p,r, provide a more sensitive search strategy for 
finding new physics in rare B decays than for example the decay B — > X s 7 , which 
constrains the magnitude of Cy . The sign of C^f , which depends on the underlying 
physics, is not determined by the measurement of BiB — ► X s + 7). This sign, which in 
our convention is negative in the SM, is in general model dependent. It is known (see for 
example ||111| ) that in supersymmetric (SUSY) models, both the negative and positive 



signs are allowed as one scans over the allowed SUSY parameter space. We recall that 
for low dilepton masses, the differential decay rate for B — ► X s £ + £~ is dominated by 
the contribution of the virtual photon to the charged lepton pair, which in turn depends 
on the effective Wilson coefficient C7 . However, as is well known, the B —>■ X s £ + £~ 
amplitude in the standard model has two additional terms, arising from the two FCNC 
four- Fermi operators |, which are not constrained by the B — > X s + 7 data. Calling their 
coefficients Cq and C±o, it has been argued in [ |111|| that the signs and magnitudes of 
all three coefficients Cj , C$ and Ciq can, in principle, be determined from the decays 
B -> X s + 7 and B -> X s £ + £~ . 

The SM-based rates for the decay b s£ + £~ , calculated in the free quark decay ap- 
proximation, have been known in the LO approximation for some time [J72"| . The required 
NLO calculation is in the meanwhile available, which reduces the scheme-dependence of 

■l- This also holds for a large class of models such as MSSM and the two-Higgs doublet models but 
not for all SM-extensions. In LR symmetric models, for example, there are additional FCNC four-Fermi 



operators involved [H2| 
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the LO effects in these decays |fL13|| . In addition, long-distance (LD) effects, which are 
expected to be very important in the decay B — > X s £ + £~ \\11A \, have also been estimated 
from data on the assumption that they arise dominantly due to the charmonium reso- 
nances J /if} and ip' and higher resonances through the decay chains B — > X s J/if>(ip' , ...)—> 
X s £ + £~ . Likewise, the leading (l/m^ 2 ) power corrections to the partonic decay rate and 
the dilepton invariant mass distribution have been calculated with the help of the op- 
erator product expansion in the effective heavy quark theory [[33 1 
have, however, not been confirmed in a recent calculation |!6] 
corrections in the branching ratio BiB — > X s £ + £~) are small. 



The results of []33|1 
which finds that the power 
Moreover, the end-point 



dilepton invariant mass spectrum in this order is not calculable in the heavy quark ex- 
pansion, which requires either resummation in the context of the HQE approach or a 
non-perturbative model. We review the salient features of the decay B — > X s £ + £~ here. 

The amplitude for B — > X s £ + £~ is calculated in the effective theory approach, which 
we have discussed earlier, by extending the operator basis of the effective Hamiltonian 
introduced in Eq. fl9"H|): 

AG 



H eff (b -f s + T ,b-> s + £ + £~) = H eff (b 
where the two additional operators are: 



s + 7) - —^-KVfb [C 9 9 + dodo] , (159) 
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The analytic expressions for C${mw) and Ci (mw) can be seen in 



(160) 



and will not 



be given here. We recall that the coefficient Cg in LO is scheme-dependent. However, this 
is compensated by an additional scheme- dependent part in the (one loop) matrix element 
of Og. We call the sum Cg^ , which is scheme-independent and enters in the physical 
decay amplitude given below, 
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(162) 



where s = q 2 /m%- The function Y(s) is the one-loop matrix element of Og and can be 



seen in literature [T(| |113|| . The dilepton invariant mass distribution in B 
(ignoring lepton masses) is, 
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with u(s) = J[s — (1 + rh s ) 2 ] [s — (1 — m s ) 2 ], and f(m c ) has been given earlier; likewise 
tv[m c ) = 1 — 2a s (/x)/37r [{ti 2 — 31/4)(1 — m c ) 2 + 3/2] is the same as the corresponding 
function in B — > X r £ve, and the rest of the functions are defined as 



a x \s,m a ) 
a 2 {s,m s ) 
a 3 (s,m s ) 



-2s l + §(1 + mi) + (1 -m 
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(164) 
(165) 
(166) 



Here rhi = mi/m b and R(C^) represents the real part of Cj . A useful quantity is the 
differential FB asymmetry in the c.m.s. of the dilepton defined in refs. ||115|| : 



dA{s) 
ds 
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(167) 



where z = cos 6, with 9 being the angle between the £ + direction and the 6-quark direction 
in this system, which can be expressed as: 



dA(s) 
ds 



3a' 



1 
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u 2 (s)C 10 \m(Cf {§)) + 2C e 7 ff (l + m 
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The Wilson coefficients Cy , 



Cg" and Cio appearing in the above equations can be 



determined from data by solving the partial branching ratio B(As) and partial FB asym- 
metry A(As), where As defines an interval in the dilepton invariant mass ||111|| . A third 



quantity, called energy asymmetry, proposed by Cho, Misiak and Wyler ||116| , defined as 

N(E e - > Ei+) — N(E e + > E e -) 



A 



N(E e - > E e +) + N{E e + > E t 



(169) 



where N{E^- > E^+) denotes the number of lepton pairs where £ + is more energetic than 
l~ in the 5-rest frame, is directly proportional to the FB asymmetry discussed above. 
The relation is 1061: 

(170) 



/ A(s) = BxA. 



Taking into account the spread in the values of the input parameters, /x, A, m t , and 
Bsl discussed in the previous section in the context of B(B — ► A s +7), following branching 
ratios for the SD-piece (i.e., from the intermediate top quark contribution only) have been 



estimated in |3q| : 



B[B 
B{B - 
B(B ■ 



> X s e + e 
■ X x r + r- 



(8.4 ±2.3) x 10~ 6 , 
(5.7 ± 1.2) x 10~ 6 , 
(2.6 ±0.5) x 10" 7 . 



;i71) 



A good fraction of this uncertainty is contributed by the assumed (±9 GeV) uncertainty 
on m t , reflecting this sensitivity as first pointed out in |[L17|| . The present experimental 
limit for the inclusive branching ratio in B — > X s £ + £~ is actually still the one set by the 
UAl collaboration some time ago ||118f , namely BiB — ► A s /i + /i~) > 5.0 x 10 -5 . As far as 
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we know, there are no interesting limits on the other two modes, involving X s e + e and 
X s r+r-. 

The leading power corrections in l/ra& in the Dalitz distribution for the decay B — > 
X s £ + £~, and the resulting dilepton invariant mass distribution and FB-asymmetry have 



been calculated in [36] and discussed in detail including comparison with the earlier cal- 



culation of the dilepton mass given in [Q. We give here the simpler expressions in the 



limit m s = 0. For the dilepton invariant mass distribution, the result is [36 



as 



-(1 - s) 2 (l + 2s) (2 + Ai) + (1 - 15s 2 + 10s 3 A 2 



a 



eff,2 



\C- 12 



+ 
+ 



1 - s) 2 (2 + s) (2 + Ai) + 4 (-6 - 3s + 5s 3 ) A 2 



effi 2 



4(1 - s) 2 (2 + Ai) + 4 -5 - 6s + 7s 2 A 2 ite(Cf A ) 



,eff v ^eff 



The (unnormalized) FB asymmetry reads as, 



& 4 - -*« 



2(1 - sfs + -(3 + 2s + 3s 2 )Ai + s (-9 - 14s + 15s 2 ) A 2 



(172) 



i?e(C 9 eff ) C 10 



+ 



4(1 - s) 2 + -(3 + 2s + 3s 2 )Ai + 2(-7 - 10s + 9s 2 ) A 2 



i?e(C 10 )C 7 eff \ .(173) 



The normalization constant Bq now includes also the power corrections in V sl-. 



B = B SL 



3 a 2 


v t * s v tb 


1 


16vr 2 




2 


f(m c )[K(rh c ) + h{rh c )/2m 2 ] 



(174) 



where the functions f(m c ) and n(m c ) have been defined earlier, and the function h(rh c 
is defined as: 



h{m c 



Ai + 



A 5 



f(rn c 



-9 + 24m 2 - 72m^ + 72m° - 15m° - 72m^ lnm c 



;i75) 



Doing the integration, one derives the (leading) power corrected branching ratio for 
B — > X s £ + £~. The decay width itself may be written in the numerical form 1|36|| : 



r HQE = r 6 (1 + ^ + c ^ ^ 



(176) 



where T b is the parton model decay width for b — > s£ + £ and the coefficients have the 
values 

C a = 0.50 and C 2 = -7.425 . 

This leads to a reduction in the power-corrected decay width by —4.1%, using Ai = —0.2 
GeV 2 and A 2 = 0.12 GeV 2 . Moreover, this reduction is mostly contributed by the A 2 - 
dependent term. We recall that the coefficient of the Ai term above is universal, i.e., it is 
the same as in the inclusive widths T(B — > X u £vi) and T(B —>■ X s + 7), but the coefficient 
of the A 2 term above is larger than the corresponding coefficient (= —9/2) in the semilep- 
tonic decay width. Hence, power corrections in T(B — > X u £vi) and Y(B — > X s £ + £~) are 
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Figure 6: Dilepton invariant mass spectrum dB(B — > X s e + e~)/ds in the parton model 
(dashed curve) and with leading power corrections calculated in the HQE approach (solid 
curve). (Figure taken from f36[| ). 



rather similar but not identical. Finally, the effect of the power corrections leads to a 
reduction of about 1.5% in the branching ratio B(B — > X s £ + £~). 

Concerning power corrections to the dilepton invariant mass distribution and the FB 
asymmetry, we would like to make the following observations: 

• The results given above in the HQE approach reproduce the parton model expres- 
sions for the dilepton invariant mass distribution and FB asymmetry in the limit 
Ai — > and A2 — > 0. 

• The power-corrected dilepton invariant mass distribution given above retains the 
characteristic 1/s behaviour following from the one-photon exchange in the parton 
model. 

• Power corrections in the dilepton mass distribution are found to be small over a 
good part of the dilepton mass s. However, these corrections become increasingly 
large and negative as one approaches § —>■ s max . Since the parton model spectrum 
falls steeply near the end-point s — > § max , this leads to the uncomfortable result 
that the power corrected dilepton mass distribution becomes negative for the high 
dilepton masses. We show in Fig. |6] this distribution in the parton model and the 
HQE approach. 

• The normalized FB asymmetry, dA(s)/ds is stable against leading order power cor- 
rections up to s < 0.6, but the corrections become increasingly large and eventually 
uncontrollable due to the unphysical behaviour of the HQE-based dilepton mass 
distribution as s approaches s mac 
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Based on these investigations, we must conclude that the HQE-based approach has a 
restrictive kinematic domain for its validity. In particular, it breaks down for the high 
dilepton invariant mass region in B — > X s £ + £~ . This behaviour is very similar to what 
has been observed earlier in the context of the photon energy spectrum in the decay 
B^X S + 1 . 

As an alternative to the heavy quark expansion, the non-perturbative effects due to the 
.B-hadron wave function on the decay distributions in B — > X s £ + £" have been evaluated 
in [ 36 1 using the familiar Fermi motion model. The dilepton invariant mass distribution 
is found to be stable against such effects over most part of the dilepton mass. However, 
the end-point spectrum is sensitive to the model parameters. The FB asymmetry turns 
out to be more sensitive to the model parameters. Further details on these points can be 
seen in 
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Concerning the LD effects in B — > X s £ + £~, it is worth noting that such contributions 
( due to the vector mesons such as J/ip and if)' as well as the continuum cc contribution 
already discussed) appear as an effective (slT^lX-^T^) interaction term only, i.e. in the 
operator Og. This implies that the LD-contributions should change Cg effectively. The 
LD-contribution from the matrix element of the four-quark operators 0\ and O2 discussed 
earlier in the context of the decay B —>■ X s + 7 can be absorbed in C7. However, as we 
have discussed earlier, to a good approximation CV is dominated by the SD-contribution. 
Finally, C±o has no LD-contribution. In accordance with this, the function Y (s) is replaced 
by, 

Y(s)^Y'(s)=Y(s)+Y Tes (s), (177) 



where lres(s) is given as [115 



Y Tes (s) = ^K(3C 1 + C 2 + 3C 3 + C 4 + 3Cr + C (i ) £ J^^"^ ' 
a Vi = M K ~ sm b ~ lM v^ 

(178) 

where k is a fudge factor, which appears due to the inadequacy of the factorization 
framework in describing data on B — > J/ipX s . With 

k (3d + C 2 + 3C 3 + C 4 + 3C 5 + C 6 ) = +0.88 , 

one reproduces (in average) the measured branching ratios for B — > J/ipX s and B — > 
i/j'X s , after the contributions from the \ c states have been subtracted. The long-distance 
effects lead to significant interference effects in the dilepton invariant mass distribution 
and the FB asymmetry in B — > X s £ + £~ shown in Figs. and |8|, respectively. This can be 
used to test the SM, as the signs of the Wilson coefficients in general are model dependent. 

It is obvious from Fig. [7| that only in the dilepton mass region far away from the res- 
onances is there a hope of extracting the Wilson coefficients governing the short- distance 
physics. The region below the J/tp resonance is well suited for that purpose as the 
dilepton invariant mass distribution here is dominated by the SD-piece. Including the 
LD-contributions, following branching ratio has been estimated for the dilepton mass 
range 2.1 GeV 2 < s < 2.9 GeV 2 in ||: 

B(B -> X,/iV) = (1.3 ± 0.3) x 10~ 6 , (179) 
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Figure 7: Dilepton invariant mass distribution in B — > X s £ + £~ in the SM including 
next-to-leading order QCD correction and LD effects. The solid curve corresponds to the 
parton model and the short-dashed and long-dashed curves correspond to including the 
Fermi motion effects. The values of the Fermi motion model are indicated in the figure. 
(Figure taken from [55]). 



with B(B — ► X s e + e~) ~ B(B — ► X s fi + fi~). The FB-asymmetry is estimated to be in the 
range 10% - 27%, as can be seen in Fig. |[ These branching ratios and the FB asymmetry 
are expected to be measured within the next several years at the forthcoming experiments. 
In the high invariant mass region, the short- distance contribution dominates. However, 
the rates are down by roughly an order of magnitude compared to the region below the 
J/^-mass. Estimates of the branching ratios are of O(10~ 7 ), which should be accessible 
at the LHC For further suggestions concerning polarization-dependent asymmetries in 
B -> XJ+£- , we refer to JTT9], pOR . 

While still on the subject of the inclusive decays B — > X s £ + £~, we note that both 
the dilepton mass distribution and the FB asymmetry are sensitive to non-SM effects. 
The case of SUSY was first studied in the classic paper on this subject by Bertolini et al. 
T2I\ . Since then, more detailed studies have been reported in the literature [111, |116| . 



For a recent update of the SUSY effects in B —>■ X s £ + £ , we refer to ||122|| in which it 



is shown that the distributions in this decay may be distorted significantly above the 
SM-related uncertainties, even if the inclusive decay rates are not significantly effected. 
Based on this and earlier studies along the same lines, it is conceivable that the FCNC 
decay B — > X s £ + £~ may open a window on new physics. This scenario is still a possibility 
even after the LEP anomaly in Z° — > bb decay has largely disappeared. The point is that 
flavour conserving and flavour changing neutral currents have very different underlying 
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Figure 8: FB asymmetry for B — ► X s £ + £~ in the SM as a function of the dimuon invariant 
mass including the next-to-leading order QCD correction and LD effects. The solid curve 
corresponds to the parton model and the short-dashed and long-dashed curves correspond 
to including the Fermi motion effects. The values of the Fermi motion model are indicated 
in the figure. (Figure taken from PEfl). 



contributions both in the SM and extensions of it. Hence, our experimental colleagues 
are well advised to measure the FCNC decay B — ► X s £ + £~ as precisely as possible. This, 
and the related exclusive decays, may turn out to be the first glance on beyond-the-SM 
landscape! 



5.3 Summary and overview of rare B decays in the SM 



The rare B decay mode B — > X s i/i/, and some of the exclusive channels associated with 



it, have comparatively larger branching ratios, 
in the SM is g ggfl : 



The estimated inclusive branching ratio 



B(B 



X s vv) 



(4.0 ± 1.0) x 10" 



(180) 



where the main uncertainty in the rates is due to the top quark mass. The scale- 
dependence, which enters indirectly through the top quark mass, has been brought un- 
der control through the NLL corrections, calculated in | 124 |. The corresponding CKM- 
suppressed decay B —>■ X^vv is related by the ratio of the CKM matrix element squared 



B(B -> X d vv) 
B(B X s uu) 



\v t 



td\ 



is 



:isr 
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\v« 



\v ud \ 
\v us \ 
\v ub \ 

\Vcd\ 

|vy 
\v cb \ 

\Vtd\ 

\v ts \ 

\Vtb\ 



Present Value 



0.9744 ± 0.0010 [0 



0.2205 ±0.0011 |J 
(3.1 ±0.8) x 10~ 3 §7| 
0.204 ±0.017 
1.01 ±0.18 " 



12 



12 



0.0393 ± 0.0028 g 
(9.2 ±3.0) x 10- 3 
0.033 ±0.009 
0.97 ±0.22 i 



10 



Table 1: Present values of the CKM matrix elements Wij\ discussed in the text and in the 



PDG review [12 



Similar relations hold for the ratios of the exclusive decay rates which depend addition- 
ally on the ratios of the form factors squared, which deviate from unity through 577(3)- 
breaking terms, in close analogy with the exclusive radiative decays discussed earlier. 
These decays are particularly attractive probes of the short-distance physics, as the long- 
distance contributions are practically absent in such decays. Hence, relations such as the 
one in ( |181| ) provide, in principle, one of the best methods for the determination of the 
CKM matrix element ratio 



V tA Vt, 83 



From the practical point of view, however, 
these decay modes are rather difficult to measure, in particular at the hadron colliders 
and probably also at the B factories. The best chances are in the Z°-decays at LEP, from 



where the present best upper limit stems [125 



B(B -> Xvv) < 7.7 x 10" 



:i82i 



Some other rare B decays such as (B° 
recently updated in |L0| and ||126 
original literature. 



B° s ) -> £ + £~ and {B°, B° s ) -> 77 have been 
to which we refer for details and references to the 



6 An Update of the CKM Matrix 

The present knowledge of the magnitude of all nine CKM matrix elements is summarized 
in Table [l]. We have discussed the experiments and theory which underlie five of the 
nine matrix elements listed there; these are the ones in which b quarks are involved. The 



remaining four are taken from the PDG review |T2j from where references to the original 
literature can also be got. Note that the value given in this table for \Vtb\ is from the 
direct CDF measurements || and not the one following from the CKM unitarity which 
gives \Vtb\ = 0.9991 ± 0.0004 [12]. Having updated the CKM matrix elements, we now 



discuss the present profile of the CKM unitarity triangle which is obtained by constraining 
the apex of this triangle given by the co-ordinates (p, rj) (see Fig. [l]) in the Wolfenstein 
parametrization H. 
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Of the four parameters, A, A, p and rj, the first two are: 



A = 0.2205 ± 0.0018 , 

A = 0.80 ±0.075. (183) 

The other two CKM parameters p and r] are constrained by the measurements of IK^/V^I, 
|e| (the CP- violating parameter in the kaon system), x d (B®-B® mixing) and (in principle) 
e'/e [AS = 1 CP- violation in the kaon system). The constraints from e'/e are not included 
due to the various experimental and theoretical uncertainties surrounding it at present 

jga. 



The experimental value of |e| is (12fl 



|e| = (2.280 ± 0.013) x 10" 3 . (184) 

Theoretically, |e| is essentially proportional to the imaginary part of the box diagram for 
K°-K° mixing and is given by ||128|| 



|e| = C\ e \B K (A 2 A 6 r/) (y c {f) ct f 3 (y c , y t ) - i) cc ] 

+ f)ttytf2(yt)A 2 X\l-p)), (185) 

where C\ e \ = G 2 F f^M K M^/(QV2n 2 AM K ), y t = mf/M^, and the functions f 2 and / 3 
can be found in Ref. [fL27| . Here, the rji are QCD correction factors, calculated at next- 
to-leading order in ||129|| (fj cc ), ||130|| (fj t t) and [|131|| ifj ct ). The theoretical uncertainty in 



the expression for |e| is in the renormalization-scale independent parameter Bk, which 
represents our ignorance of the hadronic matrix element {K°\{d^{l — 75)5) \K°). Some 
recent calculations of Bk using lattice QCD methods ||132| and the 1/N C approach [|133j| 
are: B K = 0.83 ± 0.03 (Sharpe PI ), B K = 0.86 ± 0.15 (APE Collaboration P5| ), 



B K = 0.67 ± 0.07 (JLQCD Collaboration fL36| ), B K = 0.78 ± 0.11 (Bernard and Soni 



136| ), and Bk = 0.70 ± 0.10 (Bijnens and Prades [ |133| |). The calculations given above 
are compatible with the range 

B K = 0.75 ± 0.10, (186) 
which has been used in the CKM analysis in ||. The present world average for AM^ is 



AM d = 0.464 ± 0.018 (ps)~ L . (187) 

The mass difference AM^ is calculated from the B®-B® box diagrams. Unlike the kaon 
system, where the contributions of both the c- and the t-quarks in the loop are important, 
both AMd and AM S are dominated by t-quark exchange: 

AM d = C d fj B y t f 2 (y t )\V td V tb \ 2 , (188) 

where C d = G 2 F / (Q^ 2 )M'^ r M B (f^ d B Bd ) ) \V t * d V tb \ 2 = A 2 \%1 - p) 2 + rf]. Here, fjs is the 
QCD correction. In Ref. [|130j| , this correction was analyzed including the effects of a 
heavy t-quark. It was found that f)B depends sensitively on the definition of the t-quark 
mass, and that, strictly speaking, only the product fjBiytjfiiyt) is free of this dependence. 
In the fits presented here we use the value fjB = 0.55, calculated in the MS scheme, 
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following Ref. |130|| . Consistency requires that the top quark mass be rescaled from its 



pole (mass) value of m t = 175 ± 9 GeV to the value ml(mt(pole)) in the MS scheme, 
given above. 

For the B system, the hadronic uncertainty is given by f Bd B Bd , analogous to B K in the 
kaon system, except that in this case f Bd has not been measured. The present status of 
the lattice-QCD estimates for f Bd , B Bd and related quantities for the B s meson, obtained 
in the quenched (now usually termed as the valence) approximation was summarized in 



giving 



f Bd = 170™ MeV, 
B Bd = 1.02±g : gi ±gS, (189) 

where the first error on B Bd is statistical and the second systematic, estimated by the 



UKQCD collaboration ||138|| . This compares well with a recent calculation by Gimenez 
and Martinelli, obtaining B Bd = 1.08 ± 0.06 ± 0.08 [1391. A modern estimate of f Bd B Bd 



in the QCD sum rule approach is that given in [I4T], which is stated in terms of /„■, and 
on using /„. = 132 MeV translates into 



fB d \]B Bd = 197 ±18 MeV . (190) 
The CKM fits being presented use 



f Bd ^B Bd = 200 ±40 MeV , (191) 

which is compatible with the results from both lattice-QCD and QCD sum rules for this 
quantity. The present experimental input can be summarized as 0: 



p 2 + v 2 = o.363 ± 0.073 (from \V ub /V cb \ = 0.08 ± 20%), 



(f Bd ^B Bd /l GeV)y/(l - p) 2 + r] 2 = 0.202 ± 0.017 (from AM d = 0.464 ± 0.018 (ps)" 1 ), 
B K r][0.93 + (2.08 ± 0.34)(1 - p)] = 0.79 ± 0.11 (from |e| = (2.280 ± 0.013) x 10~ 3 |192) 

The errors of the last two lines include the small experimental errors on AMd (3.9%) and 
|e| (0.6%), as well as the larger errors on mf (11%) and A 2 (14%). 

In order to find the allowed unitarity triangles, the computer program MINUIT is used 
to fit the CKM parameters A, p and r\ to the experimental values of |V^|, IV^/V^I, |e| 
and Xd- Since A is very well measured, we have fixed it to its central value given above. 
As discussed in , one can perform two types of fits: 

• Fit 1: the "experimental fit." Here, only the experimentally measured numbers are 



used as inputs to the fit with Gaussian errors; the coupling constants f Bd \]B Bd and 
Bk are given fixed values. 

Fit 2: the "combined fit." Here, both the experimental and theoretical numbers are 
used as inputs assuming Gaussian errors for the theoretical quantities. 
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f B = 2C0 + 4Q MeV , B B = 1 .0 , B^ = 0.75 ±0.10 
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Figure 9: Allowed region in p-rj space, from a simultaneous fit to both the experimental 
and theoretical quantities given in the text. The theoretical errors are treated as Gaussian 
for this fit. The solid line represents the region with \ 2 — Xmin + 6 corresponding to the 
95% CL. region. The triangle shows the best fit. (Figure taken from ||.) 




Since, there are still large theoretical uncertainties, we show here the results only for 
Fit 2. The results corresponding to Fit 1 can be seen in ||. The resulting CKM triangle 
region is shown in Fig. |]. As is clear from this figure, the allowed region is still rather 
large at present. However, present data and theory do restrict the parameters p and r) to 
lie in the range which we have given earlier in eq. (|141 ). The preferred values obtained 
from the "combined fit" are 

(p, rj) = (0.05, 0.36) , (193) 

which gives rise to an almost right-angled unitarity triangle, with the angle 7 being close 
to 90 degrees. However, as we quantify below, the allowed ranges of the CP violating 
angles a, /3, and 7 estimated at the 95% CL. are still quite large, though correlated. 



6.1 AM S (and x s ) and the Unitarity Triangle 



Mixing in the B®-B® system is quite similar to that in the B^-B® system, and the mass 
difference between the mass eigenstates AM S is given by a formula analogous to that of 
Eq. (|TB1): 

AM S = C s fj Bs y t f 2 (y t )\V*V tb \ 2 , (194) 

where C s = G 2 F / '(6ir 2 )MyyM Ba (f B B Bs ). To our accuracy \V c b\ = \Vt s \, hence one of the 
sides of the unitarity triangle, |Vt<z/AV^,|, can be obtained from the ratio of AM d and 
AM. 



AM. 



Vb s M Bs (/ 



B S B B 



AM, f, Bd M Bd (f 2 B Bd 



Vt 



ts 



Id 



(195) 



All dependence on the t-quark mass drops out, leaving the square of the ratio of CKM 
matrix elements, multiplied by a factor which reflects 5'?7(3)fl avour breaking effects. The 
only real uncertainty in this factor is the ratio of hadronic matrix elements. Whether or 
not x s can be used to help constrain the unitarity triangle will depend crucially on the the- 
oretical status of the ratio f Bs B B J f 2 B B Bd . In 0, a range £ s = {f Bs \[B~ B ~ s ) / {f Bd \[B~B~ d ) = 
(1.15 ±0.05) has been used, consistent with both earlier lattice-QCD [ 137|| and QCD sum 
rules [|141|| . Recent lattice-QCD calculations reported in [|139|| yield £;? = 1.37 ± 0.07, in 
good agreement with these values. (The SU(3)-breaking factor in AM S / AMd is Cs-) 
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The mass and lifetime of the meson have now been measured at LEP and Tevatron 



and their present values are Mb s = 5369.3 ± 2.0 MeV and t(B s ) = 1.52 ± 0.07 ps |51 
The QCD correction factor rj Bg is equal to its B d counterpart, i.e. f) Ba = 0.55. The main 
uncertainty in AM S (or, equivalently, x s ) is now f B Bb s - Using the determination of A 
given previously, and = 165 ± 9 GeV, one obtains || : 

f 2 B 

*• = {w - 6±3 - 3) (M^- (i96) 



The choice fB s \J Bb 3 — 230 MeV corresponds to the central value given by the lattice- 
QCD estimates, and with this our fits give x s ~ 20 as the preferred value in the SM. 



Allowing the coefficient to vary by ±2<r, and taking the central value for fB s \J Bb s , this 
gives H 

12.9 < x s < 26.1 , 
8.6 (ps)' 1 < AM S < 17.0 (ps)' 1 . (197) 

It is difficult to ascribe a confidence level to this range due to the dependence on the 
unknown coupling constant factor. All one can say is that the standard model predicts 
large values for AM S (and hence x s ). 

An alternative estimate of AM S (or x s ) can also be obtained by using the relation in 
Eq. (|195|) . Two quantities are required. First, we need the CKM ratio IV^/V^I- m 0? 



the allowed values for the inverse of this ratio as a function of fB d \jBB d was worked out. 
From this one gets (at 95% C.L.) 



2.94 < 



V t 



ts 



Vt 



1,1 



< 6.80 . (198) 



The second ingredient is the SU(3)-breaking factor which we take to be £ s = 1.15±0.05, 
or 1.21 < £ 2 < 1.44. The result of the CKM fit can therefore be expressed as a 95% C.L. 
range: 

V1.15/ - AM d ~ V 1 - 15 / 

Again, it is difficult to assign a true confidence level to AM S / AMd due to the dependence 
on £ s . However, the uncertainty due to the CKM matrix element ratio has now been 
reduced to a factor 5.3 due to the constraints on the unitarity triangle. The allowed 
range for the ratio AM S / AM d shows that this method is still poorer at present for the 
determination of the range for AM S , as compared to the absolute value for AM S discussed 
above, which in comparison is uncertain by a factor of 2. Both suffer from additional 



dependences on fB s y B Bs or £ s . 

The present lower bound from LEP AM S > 9.2 (ps)- 1 (95% C.L.) |J7| and the present 
world average AM d = (0.464 ± 0.018) (ps)' 1 can be used to put a bound on the ratio 
AM s /AMd, yielding AM s /AMd > 19.0. This is significantly better than the lower bound 
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Figure 10: Further constraints in p-rj space from the LEP bound AM S / AMd > 19.0. The 
bounds are presented for 3 choices of the SU(3)-breaking parameter: £ 2 = 1.21 (dotted 
line), 1.32 (dashed line) and 1.44 (solid line). In all cases, the region to the left of the 
curve is ruled out. (Figure taken from 0.) 



on this quantity from the CKM fits, using the central value for £ s . The 95% confidence 
limit on AM S / AM^ can be turned into a bound on the CKM parameter space (p, 77) by 
choosing a value for the SU(3)-breaking parameter £ 2 . We assume three representative 
values: £ 2 = 1.21, 1.32 and 1.44, and display the resulting constraints in Fig. [1^. This 
graph shows that the LEP bound now restricts the allowed p-rj region for all three values 
of £ 2 , though this restriction is weakest for the largest value of £ 2 assumed. Thus the LEP 
bound on AM S provides more stringent lower bounds on the parameters p and rj than 
those obtained from the CKM fits without this constraint: 

0.25 < 77 < 0.52, 
-0.25 < p < 0.35 . (200) 



Summarizing the discussion on x s , we note that the lattice-QCD-inspired estimate 
fs s J Bb s — 230 MeV and the CKM fit predict that x s lies between 13 and 26, with a 

central value around 20. All of these values scale as (/# s \Jb Bs /230 MeV) 2 . The present 
constraints on the CKM parameters from the bound on AM S are now competitive with 
those from fits to other data, and this will become even more pronounced with improved 
data. In particular, one expects to reach a sensitivity of x s ~ 15 (or AM S ~ 10 ps^ 1 ) 
at LEP combining all data and tagging techniques, and similarly at the SLC, CDF and 
HERA-B. Of course, an actual measurement of AM S (equivalently x s ) would be very 
helpful in further constraining the CKM parameter space. Note that the entire range for 
x s worked out here is accessible at the LHC experiments. 



7 CP Violation in the B System 

It is expected that the B system will exhibit large CP-violating effects, characterized by 
nonzero values of the angles a, f3 and 7 in the unitarity triangle (Fig. |I]) [|142j| . The 
most promising method to measure CP violation is to look for an asymmetry between 
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r(B° — > f) and T(B° — > f), where / is a CP eigenstate. If only one weak amplitude 
contributes to the decay, the CKM phases can be extracted cleanly (i.e. with no hadronic 

uncertainties). Thus, sin 2a, sin 2(3 and sin 27 can in principle be measured in Bd — > 7r + 7r~, 

{ bI — > J/ipK s and B s — > pK s , respectively. 

Penguin diagrams [|143|| will, in general, introduce some hadronic uncertainty into an 

otherwise clean measurement of the CKM phases. In the case of B^ — > J/ipKg, the 
penguins do not cause any problems, since the weak phase of the penguin is the same as 
that of the tree contribution. Thus, the CP asymmetry in this decay still measures sin 2/3. 
For Bd — > 7r + 7r~, however, although the penguin is expected to be small with respect to 
the tree diagram, it will still introduce a theoretical uncertainty into the extraction of a. 
This uncertainty can, in principle, be removed by the use of an isospin analysis ||144 |, which 



requires the measurement of the rates for B + — > 7r + 7r°, B° — > tt + tt~ and B° — > 7r°7r°, as 
well as their CP-conjugate counterparts. Thus, even in the presence of penguin diagrams, 
sin 2a can in principle be extracted from the decays B — > tttt. Still, this isospin program 
is ambitious experimentally. If it cannot be carried out, the error induced on sin 2a is 
of order \P/T\, where P (T) represents the penguin (tree) diagram. The ratio \P/T\ is 
difficult to estimate since it is dominated by hadronic physics. However, one ingredient 
is the ratio of the CKM elements of the two contributions: | V^V^/V^V^I ~ (V^/V^. 
From the fits in ||, the allowed range for the ratio of these CKM matrix elements is 



1.4 < 



V t 



id 



K 



ub 



< 4.6 , (201) 



with a central value of about 3. 

It is B s pKs which is most affected by penguins. In fact, the penguin contribution 
is probably larger in this process than the tree contribution. This decay is clearly not 
dominated by one weak (tree) amplitude, and thus cannot be used as a clean probe of the 
angle 7. Instead, two other methods have been devised, not involving CP-eigenstate final 
states. The CP asymmetry in the decay B s — > DfK^ can be used to extract sin 2 7 ||145 . 



Similarly, the CP asymmetry in B^ — >• D° CP K ± also measures sin 2 7 [146|. Here, D° CP is 



a D° or d which is identified in a CP-eigenstate mode (e.g. 7c + tt , K + K , ...). Further 
discussion on CP violation is given in ||126| , |147| , [148| , |149| , |150 |. 



The CP-violating asymmetries can be expressed straightforwardly in terms of the 
CKM parameters p and rj. The 95% C.L. constraints on p and rj found previously can be 
used to predict the ranges of sin 2a, sin 2f3 and sin 2 7 allowed in the standard model. Since 
the CP asymmetries all depend on p and rj, the ranges for sin 2a, sin 2(3 and sin 2 7 are 
correlated. That is, not all values in the ranges are allowed simultaneously. Given a value 



for fB d yBB d , the CP asymmetries are fairly constrained. However, since there is still 
considerable uncertainty in the values of the coupling constants, a more reliable profile of 
the CP asymmetries at present is given by the "combined fit" (Fit 2) p. The resulting 



correlations are shown in Figs. 11 and 12. From Fig. /refalphabeta2 one sees that the 
smallest value of sin 2(3 occurs in a small region of parameter space around sin 2a ~ 0.8- 

0.9. Excluding this small tail, one expects the CP-asymmetry in Bd — > J/^frKs to be at 
least 20% (i.e., sin 2(3 > 0.4). Note that the LEP bound AM s /AM d > 19.0 removes a part 
of the small sin 2(3 region in this tail. This is easy to understand if one recalls the relation 
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f n = 200 + 40 MeV , 
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Figure 11: Allowed region of the CP- violating quantities sin 2a and sin 2(3 resulting from 



the "combined fit" of the data for the ranges for fB d yB Bd and B K given in the text. 
(Figure taken from ||.) 
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Figure 12: Allowed values (in degrees) of the angles a and 7 resulting from the "combined 



fit" of the data for the ranges for /# d v/ Bs d and Bk given in the text. (Figure taken from 



sin 2(3 = 277(1 — p)/((l — p) 2 +r] 2 ). As seen from Fig. 0, the LEP bound removes the large 
negative-p values, which amounts to removing small sin 2(3 values. The allowed region 
in sin 2a is not affected significantly from the LEP-bound. Hence the following ranges 
for the CP-violating rate asymmetries parametrized by sin 2a, sin 2(3 and and sin 2 7 are 
determined at 95% C.L. to be 

-0.90 < sin 2a < 1.0 , 

0.40 < sin 2(3 < 0.94 , (202) 
0.34 < sin 2 7 < 1.0 . 



8 Summary and Outlook 

We have discussed some aspects of B decays in the context of SM. Flavour physics, in 
particular B physics, provides an excellent laboratory in testing calculational techniques 
in QCD, involving both perturbative and non-perturbative aspects. The applications 
presented here are by no means exhaustive but are fairly representative of the kind of 
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problems being studied in B decays and the techniques being used to tackle them. Not 
all experimental observations are calculable from first principles in QCD - this remains 
an ambitious and long-term goal. Nevertheless, the present quantitative rapport between 
experiment and theory (SM) in B decays is impressive. New experimental facilities will 
churn out a wealth of data encouraging us to ask increasingly sophisticated questions and 
seek their answers. 

A good part of B decays is accountable in QCD by virtue of the fact that the mass 
of the h quark is large enough to warrant perturbative calculations and the expansion 
parameter a s (rrib)/n) < 0.1 is small, so that leading and next-to-leading corrections should 
be sufficient. This, coupled with the working hypothesis that b quark can be treated as 
heavy, enables one to do a systematic expansion of the Green's functions in the parameter 
A/m& = O(0. 1). The resulting framework has found many applications. Illustrative of 
these are the semileptonic branching ratio Bsl{B), the electromagnetic penguin decay 
rate B(B — > X s + 7) and the average charmed hadron multiplicity in B decays (n c ), 
which are all in fair agreement with data. Some of the present theoretical dispersion in 
these quantities is expected to be considerably reduced as and when the complete NLO 
QCD corrections are available. 

The only visible question mark in inclusive B decays is the considerably shorter ob- 
served lifetime of the A& baryon, which is theoretically neither anticipated nor easy to 
accommodate. To firm up present estimates, one has to reliably calculate the mesonic 
and baryonic matrix elements of the local four-quark operators present in the effective 
Hamiltonian based on the SM. This is an ambitious calculation for lattice QCD and one 
which will probably not be carried through in this century. We have little choice but to 
sharpen other tools such as the QCD sum rules to draw definite conclusions. The apparent 
mismatch in lifetimes may owe itself to our imprecise understanding of the non-leptonic 
decays, but one can not exclude the possibility that it may after all have an experimental 
origin, like the once omnipresent (and now defunct) Z° — > bb anomaly. This remains to be 
settled in future experiments. In particular, at HERA-B and in experiments at Tevatron 
and the LHC, the A^-lifetime will be measured very precisely using fully reconstructed 
Vs. 

As emphasized in the introduction, B decays enter in five of the nine CKM matrix 
elements. The best measured of these is the matrix element \V c b\ (see Table |l|), which 
is determined with ±7% accuracy, with remarkably consistent results from the exclusive 
and inclusive decays. This can be taken as an excellent test of the parton-hadron duality 
in semileptonic decays. In exclusive decays, this precision has been made possible due 
to the theoretical developments in the context of HQET of which the decay B — > D*lvi 
remains the show-piece case. More work is needed to reach similar precision in other 
matrix elements of which two, \Vtd\ and |T4b|, are crucial in testing the CKM unitarity 
(see Fig. []]). The former, together with |Vt s |, will be measured in a variety of ways 
involving B - B° mixings and rare B decays. Present determination and theoretical 
proposals have been discussed here. Once again, the matrix elements of the four-quark 
operators play a crucial role and they have to be determined as accurately as possible. 
Fortunately, experiments will be able to put direct and model-independent bounds on 
some of these matrix elements. The case in point is the radiative decays B — > p7 and 
B — > u>j, where data on charged B ± and neutral B° decays can be used to disentangle 
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the contributions of the four-quark operators and the electromagnetic penguin operator 
with the help of isospin symmetry. Apart from testing the CKM unitarity, rare B decays 
are sensitive to new physics. The case in point here is the decays B — > X s £ + £~ and the 
related exclusive modes. Invariant dilepton mass and FB asymmetry in these decays, 
measured precisely, may reveal deviations from the SM. Such deviations, for example, are 
anticipated in SUSY models. 

Finally, the overriding interest in B decays is that they will test the CKM paradigm 
for CP violation. Present estimates of the CP- violating asymmetries predict a large value 
for sin 2/5. Since this asymmetry is measurable in a large number of experimental facilities 
being built, and there are no theoretical uncertainties in the interpretation of data, there 
is good reason to be optimistic that soon one would have first observations of CP violation 
in B decays which one can also transcribe in terms of the underlying CKM parameters, 
in particular rj and p. However, to quantitatively test the CKM paradigm one needs 
the measurement of at least one more CP asymmetries, related to the angles a and/or 
7. Some estimates of these asymmetries, related problems and possible resolutions are 
discussed in these lectures and elsewhere. The different ways of testing the CKM unitarity 
through CP asymmetries, rare decays and mixing will surely lead to an overdetermination 
of the CKM parameters, which is the goal of B physics. 

Acknowledgements: I would like to thank Hrachia Asatrian, Vladimir Braun, Christoph 
Greub, Tak Morozumi, and Matthias Neubert for helpful discussions. Matthias Neubert 
kindly provided Figure 2 based on the work in ref. |p0| - The warm hospitality of Marek 
Jezabek and the organizers of the Zakopane school is gratefully acknowledged. 

Post scriptum: 

As this manuscript was being completed, I heard the sad news of the passing away 
of Professor Abdus Salam, one of the principal architects of the standard model and 
uncontestedly the staunchest supporter of the third world science. The scientific world 
is poorer without him. For me personally he was a role model - an ideal teacher, a 
great scientific leader and a compassionate human being - bubbling with ideas, always 
enthusiastic, full of passion and free of prejudices. Alas, he is no more! These lectures are 
dedicated in gratitude to him. 
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